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INTRODUCTION 


The  US  Army  Missile  Command  (MICOM)  has  developed  an  inhouse 
capability  to  design,  optimize,  and  analyze  the  image  from  complex 
optical  systems.  A  lens  design  program  which  was  written  for  the 
National  Aeronautics  and  Space  Administration  (NASA)  was  use^d  as  the 
basic  design  program  and  modifications  were  performed  on  it  to  make 
it  compatible  with  the  SC  4020  digital  computer  plotter  routines. 

The  end  product  of  this  program  is  a  compilation  of  ray  coordinates 
and  diiection  dosines  on  the  final  optical  surface.  These  parameters 
are  punched  out  on  data  cards  and  used  as  an, input  to  a  program  which 
was  .developed  inhouse  to  do  a  complete  analysis  of  the  image.  The 
analysis  includes: 

1)  Plots  of  the  spot  diagram  and  ray  path  projections  on 
two  planes  perpendicular  to  the  image  plane. 

2)  Plots  of  the  modulation  transfer  function  of- the  system 
dc  ived  from  the  spot  diagram. 

3)  Plots  of  the  image  energy  distribution  derived  by 
scanning  an  aperture  across  the  spot  diagram. 

Several  advantages  are  realized  by  performing  ,the  design/raytrace 
operation  separately  the  image  analysis  operation  instead  of 

incorporating  the  Image  Analysis  Program  into  FOLDP  as  a  subroutine  to 
be  run  concurrently: 

1)  The  rays  from  several  runs  can  be  put  together  to  form 
a  composite  spot  diagram.  This  is  especially  helpful  in  the  case  of 
multiple  colors  or  segmented  optical  elements.  If  a  system  is  lay- 
traced  for  a  segmented  lens  over  a  wide  spectral  region*  then  the 
composite  spot  diagram  must  be  defined  by  one  run  for  each  combination 
of  color  and  lens-  segment.  The  outputs  for  all  of  the  cases  can  be 
combined  by  the  Image  Analysis  Program  to  frrro  a  spe-t -diagram  repre¬ 
senting  the  total  svs'tem  performance. 

2)  The  control  par.-meters  required  by  the  Image  Analysis 
Program  can  be  chosen  to  yield  the  most  meaningful  analysis  of  the 
final  spot  diagram.  Parameters  such  as  modulation  frequency  limits 
and  aperture  scan  step  size  can  bes::  be  determine  .  by  knowledge  of  the 
spot  diagram.  This  is  available  as  the  end  prode:t  of  FOLDP  and  can 
be  med  to  choose  the  input  parameters  for  the  It  age  Analysis  Program. 

3)  .  The  image  plane  may  be  moved  to  various  positions  relative 
to  the  final  optical  clement  without  the  necessity  of  a  rerun  of  the 
design/raytrace  operation. 


Preceding  page  blank 
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First,  the  modifications  made  to  the  Fortran  Optical  Lens  ne?jon 
Program  POIBP)  aro  described,  end.  then  the  image  l^sL  Prowls 

bili^v  of*thfleJTfG-CaSe  ^  iS  included  tc  demonstrate  the  cap a- 

1  ,1®ns  desrgn  program  .and  to  show  the  additional  irformation 

Program.8  V  Pr°CeSsing  the  sPot  through  the  Image  Analysis 

A  complete  Fortran  listing  of  the  Image  Analysis  Program  is 
included  as  Appenaix  A.  A  brief  outline  of  the  Jquired  input  cards 
for  the  Image  Analysis  Program  is  included  as  Appendix  B  t/give  the 

a  U Vn.  ^  ^  the  Pr°8rS"  3  **«*  »*«•»•  when^preparing 
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FORTRAN  OPTICAL  LENS  DESIGN  PROGRAM  (FOLDP) 

The  KOLDP^-  routines  ware  written  for  an  IBM  7094  machine  and  have 
been  adapted  to  run  on  MICOM’s  CDC  6600.  It  was  necessary  to  modify  the 
lens  profile  plot  routines  to  make  them  compatible  with  the  SC  4020 
plotter  which  is  available  to  be  used  in  conjunction  with  the  CDC  6600. 
These -routines  produce  a  scaled  drawing  of  the  optical  system  on  the 
SC  4020.  An  addition  to  the  program  was  required  to  extract  the  ray- 
trace  data  in  the  correct  sequence  of  the  execution  to  provide  inputs 
for  the  image  analysis  software.  The  modifications  an!  plot  routines 
will  not  be  shown  here  because  the  program  is  extensive  (7000  cards)  and 
the  modifications  would  not  be  applicable  to  another  computer  or  to  the 
same  computer  with  a  different  plotter.  However,  the  output  deck  for 
this  program  is  listed  and  discussed  since  almost  any  general  raytrace 
prograc.  can  be  modified  to  yield  an  identical  output,  and  this  deck 
serves  as  the  input  to  the  Image  Analysis  Program  which  can  be  adapted 
in  part  or  in  its  entirety  to  other  machines  and  plotters  with  a  minimum 
effort. 

A.  Output  Deck 

The  output  deck  consists  of  the  following: 

1)  A  card  containing  an  integer  quantity  (N)  equal  to  the 
number  of  ray.*  on  the  final  output  surface.  Format  (17.0). 

2)  A  series  of  N  cards  each  containing  the  X,Y,Z  coordinates 
of  the  rays-or.  the  final  surface  and  the  direction  cosines  QX,QY,QZ  of 
the  rays...  Format  (6E13.7). 


Some  other  data  cards  must  be  added  to  this  deck  before  it  com¬ 


prises  the  complete  input  for  the 

The  coordinate  reference 
system  convention  employed  for 
these  programs  is  shown  in 
Figure  1,  The  Z-axis  is  the 
optical  axis. 


Figure  1.  Coordinate 
reference  system. 


image  analysis  routines. 
Y 


1  FOLDP  -F  ortraii  Optical  Len»  Design  Program,  available  from  COSMIC.  Barrow  Hall, 
Univrmty  of  Georgia.  Athens.  Georgia  30601  (404-542-3205)  1  October  1 968. 
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A  ray  is  depicted  in  figure  2  as  it  appears  in  the  Y,Z  plane  for 
an  arbitrary- final  optical  surface. 
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Figure  2.  Ray  on  final  surface. 


D.  Capability 

The  FOLDP  measures  the  performance  of  an  arbitrary  optical  system 
by  a  merit  function  which  is  essentially  the  sum  of  the  squares  of  a 
set  of  weighted -aberrations.  Designing  is  a  lease  squares  process  for 
iteratively  reducing  the  magnitude  of  the  merit  function  by  linearly 
adjusting  system  parameters.  FOLDP  can  handle  systems  containing  up 
to  100  surfaces  (including  aspherics),  7  object  points,  6  colors,  and 
200  rays. 

C.  Example  Case 

A  Cooke  triplet  lens  system  was  selected  for  the  example  case. 
Figure  3  shows  the  system  before  the  optimization  is  performed. 

Figure.  4  shows  the  SC  4020  representation  of  the  system  after  the 
design-  iterations  have  been  completed. 


Figure  3.  Cooke  triplet  before  design  (off-axis  object), 
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Figure  4.  Cooke  triplet  after  design  (off-axis  object). 


The  output  deck  which  is  punched  by  FOLDP  is  listed  in  Table  !,•  and  the 
spot  diagrams  which  represent  the  intersection  of  these  rays  with  the 
image  plane  are  reproduced  in  Figures  5  and  6.  These  diagrams  are 
produced  on  the  line  printer  by  FOLDP.  Figure  7  is  simply  a  representa¬ 
tion  of  the  ray  pattern  on  the  entrance  pupil.  The  design  was  performed;, 
for  two  object  positions  simultaneously:  one  was  on-axis;  the  other 
was  1°  off-axis.  This  insures  good  resolution  over  the  intended  field 
of  view  of  the  lens  system.  The  ray trace  data  were  punched  by  FOLDP  for 
both  object  positions^  sc  the  Image  Analysis  Program  will  analyze  two 
separate  cases. 


9 


TABLE  1.  FOLDP  OUTPUT  DECK 


76 


,15218086*00-. 13696276*01-. 1876969£*00- 
*.  15218086*00-. 13696276*01-.  187t>'.'69£  »00 
.1513267E*00-. 10592876*0 i~. 112331sE*00< 
15132676*00-.  I059287e*0l-.n233l46*00 
.15072O86*OO*.7S36O39£*0O-.S763626E-al- 
-* 15072086*00-. 7S36O39£*O0-.5763626£-Ol 
.13O3318€*Og-,4SO9955£*O0-«2197705£rOl- 
-.1S03318£*OO-.4S0955SE*OO -.21977056-01 
.15014166*00-.  I5O1416£*O0-.43766t>lE-O2- 
15014166*00-.  15Q1,916£*Q0-.43768Q1E-O2 
.15014166*00  .1S014166*00-.4376801E-02- 
-.1S01416E*00  .1S01416E*OO-.43768O 16-02 
.15033186*09  .4S099S5E*06-.2)97705£-01« 
-.15033186*00  .45099556* 09-.21977056-01 
.15072086*00  .7S36e39£*O0-.5763626£-Ul- 
-«i5072C8£*OO  .75360396*00-.5763626E-0l 
.151326.2*00  .10592876*01-. 11233146*00- 
-.15132672*00  .10592876*01-. 11-33146*00 
.1521808E*00  .1369627E*0i-.1876969E*00- 
-.1521808E*00  .1369S27E*01t.1876969E*00 
.45S8046E*00-,1215479£*01-.1662033£*00' 
-»4SS8046E*00-.1215479E*01-.1662033E*00 
.45359436*00-.907iS66E*00-.lC06l27E*00- 
-.45359436*004, 9o71886E*00-.1008127E*00 
•4520884E*00~. 602784 5£*06-.5538921t-0i> 
-.45208e4e*00-i6027845E*OOri553«921E-01 
.4512U6E*00-.3008078e*00-.286i604E-ftl- 
-<4Si2ll6e*0O-.3OO8O78£*o6-.266i6~04€-0l 
.45092276*000.  -.19768801-01- 

-.45092372*000.  -.19768806-01 

.45121i6e>90  ,3008078E*OW-.286i604E-01> 
*.4S12116E*00  •  3007  <,786*00- ,28616046-01 
.45209846*00  .602/845£*C0-.55389216-01- 
-.45208842*00  .6027845£*00-.5538921E-01 
,45359632*00  .9071SS6E*09-.1008I27E.OO< 
453594 3£*00  .9071886e*30-.1008127£*0C 
.4558046E*00  .12i54?VE«ui-.1662033E*00 
-.45580466*00  .1215479E*O1-,1662033£*9O 
,7598099E*00-.1*63T34E.01-.168S797E*00 
-.7S98099E*00-.1063734E*01-.1685797E.OO 
.75663336*00-,7566333E*90-.li232l4£.30 
-,7S66333E*00-.75663336*06-. 11233146.00 
.754S982E*0&>  .4527589E*00-. 75702982-01 
-.7S45962E*00-.4527559E*00-.757029«E-01 
•7S36039E*00-.1507208£*00-,5763626E-Ol 
-.7536039E*00-.1507208E*00-.5?63626E-tfl 
.75360396*00  .1S072082*00-.S763626£-01 
-.75360396*00  .1507208E*00-.5763626£-01 
.754S982E«00  .45275891*00-. 75702982-31 
-.75459822*00  .452?569£*00-.75?02986»0i 
75663332*00  *75663332*00-, 11 233142*00 
-,7S66333E*00  .7566333E*00-.1123314E»00 
.7S9S699E400  ♦1063734£*01-,1&S57?7E*0« 
-,7598099E*0S  .10637342*01-. 168579/72*00 
.1Q6S847E*01-.9135832E*00-.19441's3E*00 
-« 10658476 »0i-.9l35832£*C0-» 19492536.00 
—106204l£*yl*.6068804£*»0-.  1477:9642*00 
-.I062341E*01-i6068804£*00-.l472964E*00 
.1059831E*01-.3028088£*00-.H92749€*00 
-•  10598316*0  l-.3028088£*00-.l  7,927496*00 
•10S9106E*01-*5375286£-l4-. 11002232 *00 


■.9979185E-62  .898126&E-01 
•9979i85£-02  .89812862-01 
.99933536-62  <699S3%7E-01 
.99933536-02  .6995347E-01 


.99930632-02 

.99930036-62 


•4996S01E-01 

.4996S01E-01 


,9911850 16-02  .29965506-01 
.96po50|E.t02  .299655OE-01 
.99851296-02  ,9985129Er02 
.99851296-62  .95851292-02 
.,99851292-02-.998S129E-Q2 
.99851Z9e-62-.9985129£-02 
•998850 IE-62 *.29965501-01 
.99B8S*lfc-62-.2996SSO2-01 
•9993003E-62-.4996S0 12-01 
.99930036-02-. 4996501 6-01 
,9993353E-d2-.6995347E-0i 
.99933S32-02-.699S347E-01 
.99791856-02-, 89872466-01 
. 9979 185£~02-.89e 12666-01 
.29955282-61  .79480756-01 
j'29955286-0.1  .79S8075C-01 
.29951932-01  - , 59963862-0 1 
129981936-01  .59963866-01 
.29478446-01  .39971266-01 
.29978446-01  .39971266-01 
.29968746-01  .1997916E-91 
.29968746-0)  .19979166-01 
.29964366-619;  1 
.29964362-010.  ... 

•2996874Er-3i-*  199791&C-01 
.29968746-6:-. 19979166-01 
.29973446-01-. 39971262-01 
.2997844£-6l-.3997i^€-01 
.29981936-01 -.59963866-0) 
•29981936-01-. 599638 66-21 
•2995528£-Ol-.T968075£401 
.299S5286-01-.7988075C-01 


•»4992258£-81 
*49922586-0 i 
*49966776-01 
.49966772-01 
•*49970156-01 
,.49970X56-61 
•*49965912-01 
.49965016-61 


.69891616-0.1 

.69691616-01 

.49966776-01 

.49966776-01 

.29962096-01 

.29982096-01 

*99930036-02 

.59930036-02 


•*499656 16-81-.99930036-02 
*49965016-01— .99930036-02 
•*4997015E-Or-,29982C9£-01 
.499 7 0 156-0 1-. 29982096-0 1 
•»4996877e-*l-.4996677£-fl 
•49966776-81-, 49966776-01 
■,*9922586-0l-.69891616-01 
*4992258£-Ol-.69691616-01 
.69837926-01  .59861086-0* 


.69837926-01 
.69923206-01 
.69923206-01 
,69949526-01 
•699495 2E-01 

-.69954566-01 


.59861086-01 

,39954116-01 

.39956116-01 

.19985586-01 

.19985586-01 

•35504086-15 


.99590876*00 

.99590876*00 

,99750026*00 

•99750026*00 

•99870106*00. 

.99870106*00 

.99950106*00 

*99950106*60 

.99990036*00 

,99990036*00 

•99990036*00 

.99990036*00 

.99950106*00 

.99950106*00 

.99870106*00 

.99870106,03 

.99750026*00 

*99750026*00 

.99590876*00 

.99590876*00 

.99335426*00 

.99635426*60 

•99.775026* 0 0 

.99775026*00 

.9967510E*00 

.99875106*00 

.99935316*06 

.99935116*00 

.99955106*00  ■ 

.99955106*00 

.99935116*00 

.99935116*00 

•99875106*00 

.99875106*00 

.99775026*00 

.99775026*90 

.99635426*60 

.99635426*00 

.99630466*00 

,99630486*06 

.99759026*00 

.99750026*00 

.99830066*00 

.99830066*00 

.99670106*00 

.99870106*00 

•99670106*00 

.99870106*00 

.99839066*00 

.99830066*00 

.99750026*00 

.99750026*00 

.99630466*00 

.99630466*00 

.99576076*00 

.99576076*09 

.99675196*00 

.99675196*09 

.99735036*00 

.99735036*00 

.99755026*00 


10 


TABLE  1.  c CONTINUED 


-,19591Q6E*01- 
,1059831£*0l 
-,10S9831E*01 
.1062041E*0l 
-.10620416*01 
.10658476*01 
-.10658476*01 
,J.37)631£*01- 
-.1371631£*01- 
,1369627E*01- 
-.13696276*01 
„1 36962 7E»Q1 
136962 7E+ 01 
< 1371631E*01 
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Figure  5.  Spoc  diagram  for  on- axis  object. 


IMAGE  ANALYSIS  PROGRAM 


The  Image  Analysis  Program  was  developed  to  perform  a  complete 
ipiage  analysis  for  any  optical  system  which  can  be  raytraced  by  FOLDP. 
Since  the  rays  are  defined  on  the  final  optical  surface,  the  user  is 
free  to  specify  the  location  of  the  image  plane  relative  to.  this  surface 
and  examine  the  image  at  that  plane.  Th  .  image  plane  location  may  be 
moved  and  the  results  examined,  for  the  ,iew  location  without  a  new  ray- 
trace  being  performed. 

A.  Main  Program 

Lhe^main  program  uses  the  deck  punched  by  the  FOLDP  routines  plus 
two  additional  cards  fisting  the  run  title  and  the  image  plane  location 
relative  to  the  final  optical  surface.  The  deck  consists  of  the 
following: 

,;/• 

1)  A  card  listing  the  run  title  in  a  field  60  characters 

wide. 

IJ 

2)  A  card  listing  the  axial  distance  of  the  image  plane  from 
the  vertex  of  the  final  optical  surface  (DELD).  Format  (F10.0). 

3)  A  card  listing  the -total  number  of  rays  on  the  final  „ 

optical  surface  (H).  Format  (HO).  (Punched  by  FOLDP). 

4)  N  cards  listing  the  X,Y,Z  coordinates  of  e*ch  ray  and 
the  associated  direction  cosines  QX,QY,QL  one  ray  per  card.  Format 
(6E13.7).  (Punched  by  FOLDP). 

The  main  program  uses  these  inputs  to  compute  the  ray  intersect 
coordinates  on  the  image  plane.  It  then  calculates  the  average  X  and 
average  Y  ray  coordinates  (centroid)  and  transforms  all  rays  relative 
to  that  point.  These  values  are  then  stored  in  COMMON  along  with  the 

centroid  coordinates,  the  direction  cosines,  the  image  plane  location, 

and  the  run  title  for  use  by  the  program  subroutines. 


B .  Subroutine  RAPLOT 


1.  Purpose 

The  purpose  of  subroutine  RAPLOT  is  to  plot  the  ray  intersec¬ 
tion  pattern  on  the  image  plane,  and  also  to  plot  the  convergence  of 
the  ray  bundle  as  viewed  along  the  two  coordinate  axes  perpendicular 
to  the  image  plane.  These  plots,  when  studied  together,  give  the 
designer  a  good  "feel"  far  how  the  image  is  formed  and  whether  the 
image  plane  is  located  at  the  best  focus  position. 
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The  RAPLOT  subroutine  uses  the  X,Y,Z  ray  coordinates  on  the 
image  plane  and  the  associated  direction  cosines  "r  construct  three 
diagrams: 


1)  The  spot  diagram 
which  represents  the  intersection 
of  the  rays  with  the  image  plane 
(X,Y  plane)  (Figure  8), 


Figure  8,  Spot  diagram. 


V 

4 


* 


-*►  x 


2)  Projection  of  the 
ray  paths  on  the  Y,Z  plane  as 
reviewed  along  the  X-axis 
(Figure  9).  The  intersection 
point  is  determined  from  the’ 

Y  coordinate  and  the  slope  is 
determined  from  the  direction 
cosines  QY  and  QZ. 


Figure  9.  Ray  projection 
on  Y,Z  plane. 
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.  3)  Projection  of  the 
ray  paths  on  the  X,Z  plane  as 
viewed  along  the  Y-axis 
(Figure  10).  The  intersection 
point  is  determined  from  the 
X  coordinate  and  the  slope  is 
determined  from  the  direction 
cosines  QX  and  QZ. 


Figure  10.  Ray  projection 
on  X,Z  plane. 


X 
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INTERSECTION  WITH  IMAGE  PLANE 
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3.  Input  and  Output  Description 
a.  Input 

The  X,Y,Z  ray  coordinates  relative  to  the  pattern  centroid 
and  the  direction  cosines  of  each  ray  are  available  to  this  subroutine 
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through  the  COMMON  storage  block.  The 
input  required  for  this  subroutine  is 
a  card  containing  the  values  for 
SCALE,  DELTXY,  NEMPHS and  NCHRS.  The 
format  is  (2F10.0,  2110).  SCALE  speci¬ 
fies  the  plot  scale  such  that  the  rays 
are  plotted  on  a  field  centered  about 
the  point  (0,0)  and  defined  by 
Figure  11.  DELTXY  specifies  L'qe  plot 
tic  increment.  NEMPHS  specifies  the 
tics  which  are  to  be  labeled  and 
emphasized.  Every  NEMPHS—  tic  is 
emphasised  and  labeled.  NCHRS  speci¬ 
fies  the  number  of  characters 
including  the  decimal  to  be  written 
for  each  label. 


Figure  11.  Definition  of  SCALE. 


b.  Output 


This  subroutine  provides  a  hard  copy  of  the  image  plane 
location,  the  ray  pattern  centroid,  and  the  coordinates  of  each  ray  on 
the  image  plane.  In  addition  the  three  plots  previously  described  are 
drawn  by  the  SC  4020  plotter. 


C.  Subroutine  MTF 


1.  Purpose 

The  purpose  of  subroutine  MTF  is  to  compute  the  modulation 
transfer  function  (MTF)  of  the  optical  system  under  consideration  and 
to  plot  this  function. 

2.  Theory 


The  modulation  transfer  function  provides  a  measure  of  the 
performance  of  an  optical  system.  Suppose  that  the  object  is  a  set 
of  light  and  dark  bands,  the  distribution  of  brightness  being 
1  +  cos  2ttu'X.  The  optical  system  degrades  the  image  so  that  the 
contrast  between  the  bright  and  dark  bands  is  less.  Suppose  the  image 
distribution  of  brightness  is  F(X)  =  a0  +  ai  cos  -tjvX.  The  MTF(u)  is 
defined  to  be  the  ratio  between  the  contrast  of  the  image  to  the 
contrast  of  the  object  for  an  image  with  the  spatial  frequency  u.  In 
this1 case,  the  contrast  of  the  image  is  ai/a0,  the  contrast  of  the 
object  is  1,  and  the  MTF(v)  = 

Smith^  has  shown  that  the  MTF  can  also  be  derived  from  the  spread 
function  of  a  point  object.  Suppose  that  the  image  distribution  of 
brightness  is  F(X,Y).  Define  the  line  spread  function  by  A(X)  -  Jf<x,y). 


-Smith,  Warren  J.,  Modern  Optical  Engineering.  McGran'  Hill,  New ’t  ork.  11)66,  [>J>.  300-325. 
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Then  the  MTF  is  given  by  the  expression: 


where 


MTF(u)  =  [a^(u)  +  Ag(u)J  ^ 


,  .  /A(X)  cos  (2T7vX)dX 

A*(u)  =  — - 


As(v) 


/A(X)  sin  (2rruX)dX 
/A(X)dX 


(1) 


(2) 

(3) 


and  if  the  line  spread  function  A(X)  is  asyranetrical*  a  phase  shift  <j> 
in  the  pattern  image  is  introduced,  described  by: 


Ac(u)  Asfy) 

cos  ♦  ■  "[awT  or  tan  *  *  vw  ‘  <4) 

3 .  •  Methodology 

Assume  that  a  spot  diagram  has  been  prepared  from  ray trace 
data,  such  that  the  coordinates  of  each  ray  intersection  on  the 

image  plane  are  known.  The  line  spread  function  a(X)  xs  represented 
by  the  number  of  rays  per  interval  N(X)  (Figures  12  and  13). 
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Figure  12.  Spot  diagram  divided  into  increments  of  AX 
to  count  number  of  rays  Nx  per  AX  interval. 


Figure  13.  Determining  line  spread  function  A(X). 


Since  A(X)  is,  in  practice,  never  represented  by  a  continuous  function, 
we  must  approximate  the  integrals  in  Eqs.  (2)  and  (3)  by  the  equivalent 
sums  over  the  AX  intervals. 


ZN(X)  cos  (2rruX)AX 
'■  Ac(l’)  "  2K(X)6X 

(5) 

,  ,  2N(X)  sin  (2n»X)AX 

•4s(0)  *  SNCxiax 

(6) 

We  may  factor  out 

AX  to  simplify  Eqs.  (5)  and  (6): 

A  M  -  SN(X)  cos  (2nuX) 

}  2N(X) 

,  .  ZH<X)  sin  (2t tv*) 

“  SN«) 

(7) 

But  2N(X)  equals 

the  total  number  of  rays  on  the  spot  diagram. 

Iv,T0.AL* 

In  the  numerators 

of  Eq.  (7),  the  sums  over  intervals  are  eeuivalent  to 

sums  over  rays. 

Thus,  Eq.  (7)  becomes 

nTOTAL 

J)  cos  (2ttvX*) 

- 

a ,.00  *  M - 

"TOTAL 

(8) 

l 


As  (v) 


nTOTAL 

£  sin  (2nuX£> 
1=1 _ 

ntotal 


(9) 


which  are  the  equations  that  are  programmed.  Once  these  are  computed 
the  MTF  is  found  from  Eq.  (1),  and  the  associated  phase  angle  is 
computed  from  the  second  of  Eq.  (4).  This  subroutine  also  includes 
a  feature  which  permits  one  to  scan  the  spot  pattern  from  two  preselected 
angles,  8^  and  62.  Tc  achieve  this,  the  rotation  of  ares  is  performed 
from  the  initial  coordinate  axes  as  depicted  in  Figure  14.  The  coordi¬ 
nate  Xj  is  transformed  according  to  the  equation 

•  4  Xj  =  Xt  cos  8  +  sin  6  (10) 

and  the  computation  of  the  MTF  proceeds  as  previously  discussed, 
utilizing  x|  in  Eqs.  (8)  and  (9). 


Y 


Figure  14.  Rotation  of  axes. 
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4.  Input  and  Output  Description 
a.  input 

The  coordinates  (X^jY^)  and  NxoiAL  have  previously  been 
computed  in  the  program  and  are  stored  in  COMMON.  The  only  input 
required  by  the  user  is  a  card  utilizing  a  (3110,  4F10.4)  format  on 
which  the  following  are  punched: 

1)  An  integer  which  specifies  the  plot  tics  which  are 
to  be  emphasized  and  labeled.  (NEMPHS:  110  format). 

2)  The  number  of  characters  including  the  decimal  to  be 
written  for  each  tic  label.  (NCHRS:  110  format). 
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3)  The  number  of  froquencies^desired  (NFREQ:  HO  format). 

4)  Tiie  desired  starting  frequency  (STFREQ:  F10.4  format). 

5)  The*  desired  frequency  increment  (DFREQ:  FIG. 4  format). 

■  ~s 

ij)  The  first  scan  angle . (THETA(l) :  F10.4  format). 

7)  The  second  scan  angle  (TH5TA(2):  F10.4  format). 

b.  Output 

Tills  subroutine  provides  a  hard  copy  output  c  *  KFRT.Q, 
STFREQ,  DFREQ,  9,  NEliPHS,  and  NGHRS  values  at  the  top  of  the  page; 
followed  by  a  columnar  listing  of  the  frequency  value,  the  MTF  associ¬ 
ated  with  that  frequency,  and  the  phase  angle.-  There  are  NFREQ  entries 
per  column,  with  frequency  values  running  from  STFREQ  to  STFREQ  + 

(NFREQ  -  1)  x  DFREQ.  Also  provided  is  a  plot  of  HfF  versus  frequency 
via  the  SC  4020  plotter. 

D.  Subroutine  SCANS 

1.  Purpose 

.The  purpose  of  subroutine  SCANS  is  to  scan  the  spot  diagram  at 
the  image  plane  of  the  optical  system  under  consideration,  using  either 
a  circular  or  rectangular  scanning  aperture.  The  limits  of  the  scanning 
cycle  are  sec  by  the  user.  Knowing  the  number  of  rays  contained  within 
the  aperture  at  each  step  of  a  scan,  one  tan  form  an  idea  of  the 
intensity  distribution  of  the  image. 

2 .  Theory 

It  is  frequently  desirable  to  determine  the  density  distribu¬ 
tion  of  ray  intercepts  on  the  image  plane  of  an.  optical  system.  This 
distribution  is,  of  course,  dependent  on  several  factors:  the  shape  and 
size  of  the  aperture  used,  the  step  size  used  for  the  scan,  the  limits 
chosen  for  the' scanning  process,  and  the  intrinsic  distribution  of  rays 
on  the  image  plane. 

3.  Methodology 

Assume  that  a  spot  diagram  of  N  rays  has  been  prepared  from 
raytrace  data,  such  that  the  coordinates  Xj,Y£  of  each  ray  intersection 
on  the  image  plane  are  known.  Ail  points  in  the  diagram  have  been 
transformed  to  make  the  centroid  the  new  origin  of  coordinates  in  the 
main  program.  Thus,  in  the  following,  X^Yj  will  refer  to  the  I—  point 
referenced  to  the  centroid.  However,  the  user  has  the  option  of  speci¬ 
fying  any  origin  coordinates  (XTR.yiR)  ha  desires.  Moreover,  if  he 
wishes  Co  rotate  the  sc.:n  c.vos  rough  a  preselected  angle  8  (counter¬ 
clockwise),  he  has  this  option  as  well.  For  example,  if  the  centroid 


is  used  for  the  scan  center  and  the  axes  are  rotated  through  an  angle 
of  8  =  30°,  the.  equations  would  he: 

XTR  =  0 

YTR  =0 

T1  *  Xi  -  XTR  =  Xi#  i  =  1,  N 

T2  ■  Yi  -  YTR  -  Yit  i  «  1,  K 

X'±  =  Tl  cos  30°  +  T2  sin  30°,  i  =  1,  N 

Y*  =  Tl  sin  30c  +  T2  cos  30°,  i  =  1,  h  . 

Note  that  the  translation  of  axes  is  done  first,  followed  by  the  rota¬ 
tion  of  axes  about  the  new  center.  We  also  see  that  when  8  »  0  that 
X{  =  Xi  and  Yj  -  Y^,  as  it  must. 

Suppose  that  have  completed  the  desired  translation/rotation 
process  and  now  have  our  spot  diagram  in  a  form  ready  for  scanning. 

We  must  now  use  some  aperture  to  scan  this  diagram.  .  The  shape  and  the 
dimensions  of  the  aperture  follow  this  convention  (Figure  15). 


Y 


Figure  15.  Rectangular  aperture  used  for  scanning  along  X-axis. 


If  the  aperture  is  rectangular,  SSI  is  the  horizontal  (X)  dimension  of 
aperture  and  SS2  is  the  vertical  (Y)  dimension  of  aperture.  For  a 
circular  aperture,  SSI  and  SS2  =  2R  where  R  is  the  radius  of  the  circl 
The  number  of  steps  NS  in  the  scan  is  determined  by  the  programmer  and 
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Figure  17.  Rectangular  aperture  used 
for  scanning  along  Y-axis. 


Note  that  here  the  scanning  rectangle  maintains  the  same  orientation 
with  respect  to  the  coordinate  axes  as  it  had  in  the  X-axis  scan.  The 
user  may  alter  this  orientation  by  inserting  a  ROTATE  card;  ‘which  will 
make  SS2  =  SSI  and  SSI  =  SS2;  that  is,  the  scanning  rectangle  will  be 
rotated  through  90°  (Figure  18). 


Y 

i 


SCAN  DIRECTION 


X 


SCAN  DIRECTION 


Figure  18.  Rotation  of 
rectangular  aperture, 
for  Y-axis  scan. 
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The  scan  Chen  proceeds  in  the  same  manner  presented  for  the  X-axis 
scanning  process,  utilizing  the  appropriate  Y-axis  parameters  in  pl*ra 
6 £  the  X-axis  parameters  previously  mentioned. 


4 Input  and  Output  Description 


Input 


The  coordinates  (X£,Yi)  and  N  have  -been  previously  computed 
by  the  program  and'-are  how  stored  in  COMMON.  The  data  cards  to  be 
supplied  by- the  user  are,  in  this  .order: 


1)  A  card  utilizing.  an  (110,  6P10.0,  A6)  format  on  which 
are  punched  the  values  for  NS, 'SSI,  3S2,  R,  YTR,  XTR,  TILT,  and  the 
KEY-,  in  the  order- given.  Note  that  if  KEY  -  CIRCLE  and  an  R  value  is 
specified,.,  this  will  override  any  SSI  and  SS2  valubs=  entered.  If  a 
circular  aperture  is  not-  desired,  R  .and  KEY  fields  .are  blank,  and 
valiid£-  for  SSI.  and  SS2  are  entered,  if  the  ;centroid  is  the  desired  scan 
center,  YTR"®  XTR  -  0,;  TILT-  is  the,,  angle  8,  in  degrees,  through  which 
the  coordinate  axes;  are.  to  be  rotated,  if  desired. 


2)  A  card  utilizing  a  (2F10.0)  format  on  which  are  punched 
the  values  for  Al  and  A2.  - 


3)  A"  card  utilizing  ah  (A6)  format  on  which  is  punched 
ROTATE  if  the  Y-axis  scan  is  to  use  the  rectangular  aperture  rotated 
through  9$°.  Otherwise,  a  blank  card  is  inserted. 


4)  A  card  utilizing  an  (F10.0,  2110).  format  on  which  are 
punched  the  values  for  DX,  the  X  increment  to  he  used  for  spacing  tics 
on  the  plotter  chart  axis,  N,  which  specif ies-which  tics  are  to  be 
emphasized  and  labeled,  and  NGHRS,  which  specifies  the  number  of 
characters  in  the  label. 


b.  Output 


following: 


This  subroutine  provides  o  hard  copy  output  of  the 


input  above. 


1)  An  echo  of  the  data  on  cards  1  through  3  described  in 


2)  A  columnar  listing  of  each  XC  value,  the  number  of  rays 
enclosed  by  the  aperture  with  center  at  XC,  and  the  normalized  value  of 
this  number  (that  is,.- the  intensity).  There  are  NS -of  these  values. 


3)  A  columnar  listing  of  each  YC  value,  the  number  of  rays 
enclosed  by  the  aperture  with  center  at  YC,  and  the  intensity.  There 
are  NC  of  these  values. 


.  X  tuu.*3L*  ■ 


t*'*'*  *p.-*r+u  >/„  '*  s*.  t^-r 


k«--  -  *rt**r**y~: ,  ^rf  ^/»^r«r^'>^‘V^>  *. 


4)  The  ^values  DX,  N,  and  NCHRS  specified'  in  card  4  ;.i»  input 
above.  Also  provided  is  a  plot  of  intensity -versus  XC  and  YC  values 
via  the  SC  4020  plotter. 

E.  Example  Case  Input  and  Output 


The  Cooke  triplet  raytrace  data  -which,  were  punched  out  by  TOEBP 
are  used  as  part  of  the  input  deck  for  the  Image  Analysis,  Program*  The 
complete  input  data  deck  for  this  case  is  listed  in  Table  2,  gHe  output 
from' the  line  printer  is  listed  in  Table  3,  and  the-  graphs  from  the 
SC  4020  are  listed  in  Figures  19  through  28. 


TABLE  2._  INPUT  DECK  TO  THE  IMAGE  ANALYSIS  PROGRAM 
;  FOR  THE  EXAMPLE  CASE 


08-AXIS  OBJECT  c  EXAMPLE  CASE 

•  15.  ' 

.  yb 

..,,S2I8(f9';*00<**ia69£272*01-.  18769692*00- 
15218082*00.-. 1369627F*01~, 18769692*00 
.  l5i32fi?.E*00-»  10592S7£*0t~.  11233142*00- 
15132672*00-.  1059297£*.01-»1  1233142*00 

♦  !5t)72O8E*P0-.7536O39E*O0-‘i£?63fi26!->dl- 
t ^6 ^2C8E‘ O.a^.  75360 39‘fc+ 00- . S 7636£6£-0 1 

.15O33l9E*OO-,4S09955E*0O-.2l977G5E-0i'- 
-.1503319E400-. 45099556*00-. 21977052-01 
*1-551 41 6£*00-*  !.5O1416£«OO-.43768OlE-.02- 
-.15O1416£*bO-.l5O1416E*OO~.43748OlE;<-02 
.1.E01416£*00  « 19Ol4I6E*00~»V3768QlE“02^ 
-.15014166*00  .l5wl4l6E*OQ.-*4376801E-Q2 
.150331«E*00  .4S09955£*00-.2i97705£-01- 
-.15(J33iaE*00  .450995SE.*00-.2197705E-Ql 
•1507208£*00  . 7536039E*00-, 57636262-01- 
15(/7209EtfOO  .75360392*00 '.57636262-0 l 
.m-3267£>oo  *io:r9287£*Qi-.u233i4E*Qo- 
-.I513267£*b<i  ;vl05923?E+oT>-  >11233142*00 
•  .15218086*00  .i369627E*01-i:976969£i00- 
r.lS21?08£*fi0  .136?621E*91-.l876969t*00 
»4S  58?46£*fl0-il2l5479£*0i-. 1862033E+00* 
~*4558346E*6or.m34?92*01-.  16620336*00 
».4'93S9J,.3E+00-«907tii36E+QO-\  1.0081276*00- 
r'.4535943£*bO-.907183AE*00-iT008127E*00: 

,-45208846'+00-.  60278  45E+Q0-,  552392l£-0 1- 
-.4520884?  *  00- ■i6t>27845E*0(»-.  5538921 6-01 
-  -v45T2U6E*00-.3009078C»00-.2861h04£-01- 
-*45121  l6E*0O-«3Ob3O78f>0O~,2fl6l604£-Cl 
.<♦5092  376*000.  -.1976830E-01- 

-.45032372*000.  \  ~* 1976B80E-01 

*45121 186*00  .3008a78£*,00-i286l5d4c-AI. 
-*45l2U6£*00  .300Sd?8E*G0-.28616Q4£-0l 
.4520384E*00  .602?845E.*09~*  5S38921E-01- 
-*45?9;88A2*O0  »6027845E+00*-«SS38921E-Oi 
*4S359Vit  *00  ,  8r/71886E*00~.  1008 1276  *00- 
-.45359436*00  * 907 18B6E430-. 19081 276*00 
,45,59046:1  *00  ,1215479£*014.i662033EtOO- 
-.45590462*00  .12154796*01 -.16620336*00 
. 75980496*00-. 1 •63734E*01-.168579?6*00- 
- »7'5C 9099F * 9 j- *10637342  »91-,168S797?.*00 
,7966J33E*6G-,7566333E*00-.l 1233146*00' 
-*7.566333t  *OO-.?566333E*O0-.  1 123314E  *01! 

. 75459826*00-. 45275B9£*O0-I757O?98E-Cl' 
-.754593?E*.p.v-»4527599E*O0-.. 75702982-01 
*75  560 39£*OO-*15O7208£  +  U0-.576362o£, -{.'!■ 
'  .  75  36  0  39E  *  00  -» 1 59  7"?  Ofifc  ♦  0  0  -  .  57*3626?  -  0 1 
*75360392*00  .1507208£*00-.5?63626E*31- 
-,?$36P39**00  ,  1F.07208£*00-. 57636266-01 
.75459826.00  , 452 75B9E.O0-, 75702986-01 
-.75459855 *00  *45?7589£  *00-. 75702986-0 1 
>  -  .7'jr»6331E60O  ,7S683.23£*00-.1123314E*00 
-.7066333- *00  ,.756s333f;*,00-7  1I23314£*00 
»  7«98099r  ->  ftp  *  1  66 37 34£  *  0 ;  - .  1 68579 ?E  *00 
-.75980992*09  < 1863734£*01-.1685757£ *00 
.1968847E*01-»9135832>  *00-.  19.491532*00 
-  .  X065R472  *  0 1  -  *  91  .V 58  >2£  '*00  -  .  i 949 ).  5  3E  *  0  9 
.  10620>«  'i£40i.-»6O65H045-»au...  1472969E  »08 
-.iO.6204l£*«t-.6OC3J>05»E-*OO— ;/4?2564F.*C0 


-.9979105E-82  .S981266E-01  .99590S7E*OO 
.9979185E-02  .8981266E-01  ,9959087E*C0 
.9993353E-02  .6995347£-0l  *9975002E*66 
.9993353E-02  .6995347E-01  ,9975002£*00 
■i9993Cfl3E'-02  .4996501E-01  .998701 0E*00 
. 99930.0 3E-02  .'4996501E-01  . 9.9870 10E* 00 
•»? 1885012-02  .2996550E-01  , 99950 1QE*00 
.9988501E-02  .2'996550t"01  ,9995010E*00 
-.9985129E-02  .9985129E-0?.  .9999003E+00 
i9985l29E-02  .9985129E-02  .9999003E*00 
99851292*02-. 99.85129E-02  .9999OO3E*O0 
;*9955i29£-02-,9985T29E-O2  ,9999003E*0P 
>, 998850 l£r 02- ;299C550E-O1  *9995010E*00 
. 998850  l£-0«l-i2996550E-01  .999S010£*00 
-»999.3003E*r02-«  4996501  E-OX'  .9987010E*0C 
.99930  03Er02- . 499650 1 £-0 1  .99870 1 OE+OO 
-.9993353E-Q2-,6995347E?01  .99750C2E*00 
.9993353E-02-.6995347E-01  .99?5002E*60 
•  *  9,9791 85E-fi2-,.89Bi266E-0X  *9959Q87E*00 
.997'' 1852405-. 8983 266E-01  .99590B7£*00 
-^2495528E-01  .7988075E-O1  ,9963542E*00 
.2995S28E-01  .7988075E-01  .9963542E+Q0 
-.2998193E-01  .5996386E-0P  .9977502E400 
t299H1.93£-01  .5956386E-01  *997?502E*00 
i-»?997844E-01.  .39971 25E-01  ,9987510Ei00 
*29‘»7344E-0i  .3997135E-01  .99875105*00 
-.2996874E-01  . 1997910E-01  -„99935Tl£*00 
.2996874E-01  .1997916E-01  .99935liE*09 
-.2996436E-010.  -  .5975<510E*00 

.2996436E-010,  .  .959551 OE* CO 

» *2996874£- GX-. 19979 16E-01  , 99935112 *00 
.2096874E'‘Ol-,15979i6E-Ol 
-.2997844E-C1-* 39971 26E-01 
. R997844E-0 1 -.3997 1 ?62>di 
29531 93E-01-.5996386E-01 
.  29981 93E-0 1  -.  59963.36E-0 1 
-.2995S26E-01-.798BQ75E-01 
. 29955206-01-. 7988073E-01 
-.4992258E-01  .6989161E-01 
,4952?50E-01  .6939141E-01 
-.49966772-01  ,49966772-01 
» 4  ?9o6  77E  *0 1  .4996677E-01 
-.49979)  5E-03  ,29932091.-01 
.4997015E-0}  *29982092 »01 
-.4996501 £-01  .99930O3E-O2 
*.4:99659 1E-0)  .9993003E-02,' .9581l)lOE+00 
-.4996501 E-01-. 9993O03E-02  , 998701 0E*OO 
.4996SO1E-01-.9993O03E-O2  .998701 GE*00 
», 49970 15E-01-.2998209E-01  .9983006E*00 
.  49970.1 5t*Ci-»29982O9E-01  t.9983006E*00 
-.499667VE-O1-.4996677E-01  i5975OO2E»0O 
.  499667.7E-O1*54996677E“01  .9975002E400 
••«4992268£~01-,69391£lE-01  ..9963046E400 
,49922%8E-0)-.698916lE-Ol  ,9963046E*O0 
-.63837922-01  .50851O8E-O1  .9957607E*A0 
.6983792E-01  ,59861032*01  *9957607? *00 
-.09-723202-01  .39956112-01  *996751'3E*00 
6i-99?320E-Ol  .399561 12-01  :,9O6?S19E*O0 


I 

i 

l 


,999351 12+00 

! 

.9987510E400 

_ 

.958?STOE>0O" 

-  '  |  ’* 

-  - '  M 

.9977502E+00 

'  , 

,9977502E*00 

' 

jbm 

,9903S42E*00 

1 

.9°'-3542E+00 

n 

i9'-o3046E*00 

- 

M 

f9963086E<-00 

.9975002E+00 

*99750921^00 

. 

.99,830062*00 

.99839C6E*00 

M 

,99370 102*00 
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TABLE  2.  CONTINUED 


.ios9S3is:«ai 
-.}Q598315*C1 
.  105910,6£*bl 
-.10591065*01 
« 1 05993 1E*01 
•♦10593315*01 
.10620415*01 
r. 10620415.01 
.10653475*01 
-.10659475*01 
.13716315*01 
-.13716315*01 
.13696275*01 
-.13696275*01 
.13696275*01 
-.13596275*01 
.137i'631£*01 
-.13716315*01 
.003 

4 

>-  21  . 

-  003  ~ 


. 30200595*00-. 11927495-00-. 69949525-01 
.30230995*0,0-.  11927495*00  *69949525-01. 
-.S37S296E-14-.1 1002235*00-, &995456E-01 
-.53752365-14-.  11.002235*00  ,69954S6g-01 
.30230985*00-. 11927495*00-. &9949525-01- 
.30230985*00-. 31927495*00  .6394952E-0 I* 
. 606980,45*004. 14729645*00-.  6502320E-01- 
«60698045*00r. l472964g*00  .59923205-01* 
.91 358325*00-.  19491535*00,-.  69837925*01- 
.91358325*00-. 19491535*00  .6983792E'0T- 
-,4572lb45*00-.2070071E*00-. 39752665-01 
-.45721 0^5*00- .20,70071 5* 00  .8.97S266E-01 
-.  15216085*00-.  137^9695*00-,  3981^65^01 
-.15218085*00-. 18769695*00  .39812665-01. 
.i521«08EvOC-„1876965E'*00-.999.\266E-01- 
,152180ftc*00-. I “769695*00  ,898I2$6E-0I- 
.-45721 045*90-,  207007J  5*  00-,.  39752665-0 1- 
'.4572; 045*00-. 2b.76671E*5«  ^.99752665-01- 
0002  IS  4 


.19985535-01 

,19985595-01 

•  3550.408E-/15 

.35504085-15 

.19985595-01. 

.19985585-0} 

.3995X115-01 

,39956115-01 

1.59861095-01 

•5936108E-01 

.2991755E-01 

.29917555-01 

.99791855-02 

.9979185E-02 

c9979ifl5E-b2 

*99791855-02 

■.29-#  17555-01 

•.2991 7555-01 


•9973503E+00 
.99135035*00 
.99755025*00 
.9975502E*00 
,99735035*00 
.9973,50.35600 
.9967S*9E*0C 
.99675195*00 
,995 76075 *00 
.99576075*00 
;9955146E*G0 
.99551465*00 
,99590875*00 
.99590875*00 
.99590875*00 
.99S3p8?E*00 
•995SI46E*00 
.99551465*00 


19 


50." 


.0 


90c. 


3015  : 
.003 


.00,15 


5 

EXAMPLE 


.0002 

O^r-AXlS  -04J5CT 
15.  ' 

76  ... 

.1926931E*00-.,127T356E*01- 
-,15209315*00.*.  1271 3565*014 
.  15126665*004, 96v045B55+00- 
15126665*00-. 9604,5955*00- 
. 15068295 *0b-,«54j430Eibb- 
1506829e*b0-,65434305*00r 
Vl5031335*00-.S314S19E*00- 
-.I50311?5.0b-.3514519E*00- 
. 1 50 14 11 c ♦ 0  0 - .50474 705-01- 
150141 15.00-. 50474705-rOi- 
.1501,590.“ *00  .24977035*00- 
-.1501590E*o0  .244?7d3F*0a- 
.  1'5!)36'79E>00  .550427.35*00- 
-.15036745*00  .55042735.00- 
.15077775*00  .85260905*00- 
-.15077775*00  ,eS26830E*00- 
.15} 40625*00  ,11578305*01- 
-.15140625*00  , 1 157S90£*31» 
.15229415*00  ,1.4676595*01- 
-.15229415*00  ,14<v76595*0l- 
. 45557995*00-. 1 1171475*01- 
-.45557995*00-.  11 17.1 4?6*0l- 
,453445S5*00r.8083599E*0b- 
-. 45344555*00-. 80835995*00- 
, 45200285*00-. 50359655*00- 
-.4520026E<  00^*50359565*00- 
,451 1829E+00-.201 41505*00- 
451 18295*00-, 20141 505*00- 
.45094945-00  .99434785-01 
-.45034945*00 
.45124295.00 
-.651?929E*00 
.45222955  >C0 


4, 

CASE 


♦05434785-01- 

t4001250E*00- 

,40012505*00' 

*70132995*00' 


.16162616*00- 

.161626 15 *00 

.92567236-01- 

.92567235-01 

,43938745-01- 

.4393874c -01 

,14197936-01* 

.14197986-01 

.74352566-02- 

.24352566-02 

.9 24 824 7c- 02* 

,82482375-02 
.31691066-01- 
.3i691066fOL 
.73280995-01- 
.73280995-01 
.13406815*00- 
.1 34068lE*00 
.21579065.00* 
.21574065+00 
.1432360E+00* 
,14323606*00 
•.34O4650E-01' 
■.84046506-01 
•  , 44628386-01. 
•.44&2838F-01 
■.23744956-01 
.23744955-01 
-;20734295-t)l 
■.20734295-01 
■.3542068E-01 
•.35420686-01 
'.6809249E-01 


■,99530265-02 
.99,930266-02 
■♦1 0000325-0! 
,3.000032^401 
.99962316-02 
,999.62315-02 
.99897136-02 
.598.97l3c.-02 
-.99851795-02 
.99851 79t>02 
■.93842375-02 
..  998423  7E-0? 
•.99863195-02- 
99863 1 9E-;02- 
-.99885295-02- 
.99B8529E-02* 
•.9984653E-.02- 
.99846535-02- 
99627895-02- 
.99627895-02- 
-.2998324E-01 
.29988245-01 
-.C999840E-01 
.2999860E-01 
-.29985025-01 
.2998582E-01 
-.29970956^01 
.29470956-01 
-.29963055-01 
.29963055-01 
-.29963725-01* 
,29963725-Oi- 
-.2906766E-01- 


.10074455*00. 
.10074455*00. 
.80802686-01 
.80R0269E— .01 
.6073204E401 
.60782045-01 
.407-731 95-01 
. 407731 9E-01* 
.2079277^01 
.2079277ErOl 
.82635895-03 
.62635895-03 
.1914565E-01 
.1914565E-01 
.3912549E-01 
.39125495-01 
.59067345-01 
,59067345-01 
•.7e821 996-01 
.7882199Er01 
.90774015-01 
.90774016-01 
.7079805E-01 
.7079805E-01 
.50743565-01 
.50783566-01 
.30786155-0! 
.30786155-01 
.10808725-01 
.10808725-01 
>.9164072E-02 
-.91640721-02 
-.2914142E-01 


».9948621£*00 
.99686215*00 
.99667996*00 
.99667935* 00 
..998T010E*00 
.99810105*00 
,99911855*00 
.99911355*00 
.99973395*00 
,59973395*00 
.95994985*00 
.95954985*00 
,9997fe68E*00 
.9997668E*00 
,99918445*00 
,9991844E>b0 
.99820415*00-  •> 
.9?82041E*00 
.,99683B9E*00 
,9968389E*00 
.99541995*00 
.99541995*00 
.99703955*00 
.99703955*00 
,9982594E*00 
.99825946.00 
.99907665*00 
.99907665*00 
,99949266*00 
.99949265*  00 
,999S0995*00 
.99950905*00 
,99912605*00 


29 


tU.  Sail 


TABLE  2 .  CONCLUDED 


.45222956*00 
.45380456+00 
.45330455*00 
•4S610056+00 
.45610056*00 
.75943616.00 
.75945616*00 
.75643066*00 
.75643066*00 
.75449936*00 
.75449936*00 
.75360096*00 
.75360096*00 
.75369576*00 
.75369576*00 
.75473946.00 
.75478466*00 
.75693636*00 
.7569.1636  *00 
.76024776 *00 
>7602477 f *00 
.  1 0654536"  *01 
.10454536*01 
.10613156*01 
.10519156*01 
10597546*01 
1059?c4F*m 
10591.16*01 
10591716*01 
10400406+01 
10600406*01 
10424076.01 
1 052407- *0 1 
10464'<;oE.01 
10664006*01 
13714166*01 
13714166*01 
13496176*01 
13696176*01 
13694146*01 
13699136.01 
13720356*01 
13720356*01 
.003 

4 
21 

-.003 


•701828RE+00 
.10058186+01- 
.10058136*01- 
.131 3580E+01- 
.13135906*01- 
-.96560536*00- 
-.96560536*00- 
>-.6580/ 746*00- 
.65807746*00 
.35392076*00- 
-.35392076*00- 
-.5175816E-01- 
-.51753166-01- 
.24964586*00- 
.2496469E+00- 
.55149036*00 
.551490.16*08 
.85502176*00- 
.85502j 76*00- 
.11616546*01- 
.11416546+01- 
-.81593626*00- 
-.91593626*00- 
-.50888036*00- 
-.50898086*00- 
-.20460546*00- 
-.20440546+00- 
.93248476-01- 
.99248476-01- 
.40094326*00- 
.40094326*00- 
.70474996*00- 
.70474396*00 
,1011f)53E*01- 
.1011053E+01- 
*.36006396*00- 
-.360 Ob 396  +  00- 
-.54909326-Oi- 
-.5490992E-01- 
.24941976*00- 
.24941976*00- 
.554?642E*00- 
.55426426*00- 
0002 

3 

0015  ,00i 

.003 


.68092496-01 
.11954056*00 
.11954056*00 
191 1707F+00 
.19117076*00 

•  148S9316*  00 
.14859316*00 
.98493706-01 
.93488706-01 
.67842126-01 
.67842126-01 
.55669786-01 
.55669786-01 
.61538156-01 
.61538156-01 
.85505726-01 
.85505726-0!: 
.1281317E+00' 
.12813176*00 
.19055566*00' 
.19055566+00 
A  7782786*00- 

•  1 7782786*00 
.13633796*00- 
.1 3633796*00 
.11431  ’36*00- 
.  11431 136  +  G0 
.11009326*00- 
.11099326*00 
.12613426+00- 
,126i842£*00 
.16021406*00- 
.16021406.00 
.21399036*00'- 
.2130 J03F+0O 

19892726*00- 
.19892726*00 
.1856659C.00- 
.18366596-0*5 
.19168386*00 
.19168386*00 
.21705406*00- 
.21705406*00 
S  4 

19  .0 

S 


.29967666-01 

-.29960896-01 

.29960896-01 

-.29915776-01 

^2991 5776-0 i 
-»499T1 126-01 
. 4997 1 12E-01 
-.49991 18E-01 
.4999118E-01 
-.49980586-01 
.49980586-01 
-.49966446-01 
.4996644E-01 
-.'♦9953746-01 
.40958746-0,1 
-.0995416E-01 
.4995416£-0i 
-.49935236-01 
. 4993523E-0 1 
-.49864166-01 
.49864166-01 
-.69904306-01 
.6990430E-01 
•.6995649E-01 
.69956496-01 
*59961756-01 
.69961756-01 
•.6595073E-0 1 
.59950736-01 
.69929CRE-01' 
.3992908E-01- 
.'>9879946-01- 
.69879946-01- 
.6975823E-01- 
.69758236-01- 
.39795466-01 
.39795466-01 
.39821756-01 
.39321756-01 
.89791526-01 
.39791526-01 
.89696106-01- 
.3969610E-01- 


•*•.29141426-01 
.49099666-01 
.49099666-01 
-.68939216-01 
-.63939?1E-01 
.80767236-01 
.80767236-01 
.607951  76i-01 
.60795176-Ci 
.40793176-91 
.4079317E-01 
.20800706-01 
.20800706-01 
.81892246-03 
.81892246-03 
-.29157^36-01 
-.1925763E-01 
•»  391 13S3E-01 
-.39113536-01 
-.58975056-01 
-.S397505E-01 
.7073276E-01 
• 7073276E-01 
.50785176-01 
. 50785 17E-01 
.30799136-01 
. 307991 3E-01 
.10812616-01 
•  1031 261E-01 
.91625166-02 
.91625166-02 
.2V10640E-01 
.29106406-01 
.489533X6-01 
.48953816-01 
.40755526-01 
.40755526-01 
.20803456-01 
.20803456-01 
.85039906-03 
.8503990E-03 
. 1906269E-01 
.19062696-01  , 


SO. 


.0 


90. 


.0002 


30 


-  W  1  -  T 


r  ii..  -i-ii,  - 


.99912606500 
.99834446+00 
•99834h4E*00 
.99717226*00 
.99717226+00 
.9954795E+00 
.99547956*00 
.99689706*00 
.99689766*00 
.99791686+00 
.9979168E+00 
.99853436*00 
.99853436+00 
.99075096+00 
.9987509E*00 
.9985678E+00 
.9985678E+00 
.99798636*00 
.99798636*00 
.9970,1  33E+00 
.99701336*00 
.99504286*00 
.99504286*00 
.99625656*00 
•9962565E+00 
.99707416*00 
.99707416*00 
.99749186*00 
.9974918E+00 
.99750996*00 
.9975099E+00 
.99713076*00 
.99713076*00 
.99636206*00 
.996362CE+00 
.9951260E+00 
.99512606*00 
9957406E+00 
9957406E+00 
99596026*00 
99596026*00 
99578676*00 
9957867E+00 


TABLE  3. 


LINE  PRINTER  OUTPUT  FROM  IMAGE  ANALYSIS  PROGRAM 
FOR  EXAMPLE  CASE 


0N-AXIS  39JtCT 


CX41PU  Cftse 


IHAGC  PlAN£  0IS709CF 


.  I533)33£*l? 


Y  AV£P*(,E 


•,Zht>3t 1£  05 
.26653316-05 
-.746J539E-3* 

.  44C-34 - 

.640?961£-‘4 

-.54379*16-34 

.20955646-03 

-.28*41646-33 

.JC612356-03 

-*30.>82Z«-44 - 

.J163235E-3J 
-.3{>6ft2a«-34  - 
.239*5646-03 

-.20945*4f~4* - 

.54379616-J4 

-.$4  87**  **-*4 - 

-.74635  395-04 
.74635396-04—  ~ 
-.266$3C!E-0< 

-.16 C64276-0J 
.1694437, 

-.1  7310656-03’ 
.t7A.369r-a*-- 
.U6492«£-CJ 
-.19649296-04. 

*51046626-03 
-•5 304*6 JC-. 

•65?>7  29E-CJ 
-.66  377496-94 
.5J:46S3f-?3 
-.5>C466J£-41 
.ISM92*C' JJ 

-•*7  410656-0  3 
.17M0*5C>'04- 
-.t*56427f-23 
.1656427<-8J 

-.esj47i7r-cj 

»ZS  147176-04. 


-.27345276-03 
.37  Jl^S^76-0a- 
-.:9336J!£-94 
•  l?i;t4l£-0*~.  « 
.2703734f-0.1 
-.22447446-44.... 
.27017.146-01 
-.27047346-04.  - 
-.193Jf.3lC-04 
.19336316-04 
-.37H527f-3J 


.24m<.6e-j:v 

•?*Ml4*C-J4 

.5221J2SE-03 

- ,942447*4-02 

-.27037346-03 
-,27027  3*6.43 
-.62971426-93 
-.67971426- 03 
-.39692356-03 
— .*40*02356-93 
. T3( 32  356-13 
-• 10(92 J5E-CJ 
.62971426-03 
*  -  .62971 4?6-C3 
.27337346-33 
— —  24447,346-83 
-•5&2H3266-03 
-.52244266-33 
-.2*111466-04 
.24111*66-04 
.44I69/9C-03 
♦*1*9796-03 
.15621296-03 
• .15621296-33 
-.24935676-0* 
-.2*995676-33 
-.1SJ71C5E-C3 
-33 

C.  -  * 

c. 

.35371046-03 
.153/1096-44 
;C491567£-33 
i  .2403567^-3  1.  - 
-•  15(21 29C-0T3 
-.45i.4t29C.-C3' 
^.44169796-34 
-.*41*9796-03- 
•35469006-03 
-*464*92  0£*ft3— 


•  .253471 7£-0  3 
.25347176-03 
.16793276-33 
..16743276. 53 


I7J%527E-C3 

J 744527*. OJ - 

11 621786-04 
~  --41*21 74  *-*4  - 
-.56C7961C-04 
-~~5*449*lC~4*~. 
.54079616-04 
-.S4*79*U-C4 
- • 11621 75t-C4 
•  — .11*217*6-3* 
..17345276-03 
-04 

-• 15«»*93Cr*0J 
-.J54C9006-03 
-.12674176-31 
-.12t/4i7£-CJ 


31 


.  ir4a  O/uWf 


TABLE  3,  CONTINUED 


*.5570%33£-G3 

.557*4301-31 

-.56675245-33 


.3144*52093 

.3i«4*57C-C3 

.1&1659J033 


-.5667524033  .16165*3033  3. 

.5447S*6£«*X - .161*643083  0. 

-.5047961033  .2563040017  fi. 

.53478610*3-.  .  .25633*3017  3. 

-.566?SZ6e*93  -.1614553003  C. 

—  I  1 

-.557O630C-SJ  -.»X*9*5J£-03  0. 

- — >i.4)Uli2i*0)  9, 

.lwsazrc-oj  .i^67<.srt-oj  «. 

- SJk7V17£-4J  «  J. 

.»K«<)J*f-eJ  •  .JJtll.m-tJ  c. 

--.V.1!»33*£..J3 - .  .  .0.  .. 

*?&&5 JCl£-05  C. 

'  .Z.Uiilt'Cl. - _.it4S3Slt-«t_ _ i_.. 

-.aumif-Jk  -.26t5J3lf-35  3. 

.wiitu-u - _ a.  _ 

•9l2f-(J9f-03  .  JUm'jE-CJ  0. 

» »28U.UW*fil  a  , 


PLOt  St££ 

- 


r:t  s.tetM;  nc  £*M>*»smo  ut  c«s»jci£«5 


NF££3«  1«  Sf^.ftHC  fftEO*  9,050c 


cur*  r>:a«  jc .saia  r"*«et  o*:oui:cKt 


«.  rxc  cmPACTCfi 

z*  . 

F9£0U'4C» 

ntf 

****SC  IS'^L! 

6 « 5C5 

i.ooat 

C. 053  J 

50.480 

.99  6t 

f.3<33  3 

135 .:ca 

.9776 

)•  333  3 

153.33) 

.9532 

3.0033 

.  288. C63 

359. aco 

.$665 

C.SS3 ) 

335.590 

.5125  - 

C.CCS  3 

353.333 

.7521 

3.53)3 

4QC.eC9 

-  ■  .6*63-- 

0.3030 

tVO.COS 

0.9383 

555. 353 

.<.793 

:. j;i) 

655.333 

-  ..  .*116 

C.9C53 

550.039 

.3672 

8.0033 

7S3.83C 

.2*75 

0.SC33 

756.303 

.2137 

5.999) 

333*333 

- 0666 - - 

5  •  3  3  3  Jk- 

550.300 

.1670 

C.3C33 

930. t:c 

.1156 

S.CC23 

32 


i/JKax. 


rl  Thi  i--  l*  >-  -•if  rfi  r  A-'-;  *1'-  -»  A*-  >  waiu%  . 


nc.  z?  »i  *..**«'.*  n*:**  .m»  .cm  <*i;ius*  -i.jcoc 


r  T*>6K5Uit(*  -4,?”S0  '  «36.;:ki?f*  :  Ilf  *uSW>'  -9.3)50  A»f«T<K£  <£▼* 


■  irri  Wit  »i  '■.*4n»  .;)•« 
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TABLE  3*  CONTINUED 


*  VUUC 

no.  c r 

-  f*tC<Ll 

-.0030 

3 

2.0003 

- *  -.«8?7 

9 

4.0900 

-.142*» 

> 

3.434} 

- .«2i 

•4 

O.I80C 

-.aate 

• 

3.0306 

8 

<3.8089 

»e 

•  UU 

-.3349 

?* 

.J*.2l 

-.3694 

66 

.4570 

-.3434 

46 

.468*. 

-.0303 

tH 

t.osoe 

- -  -  *W«- 

ht> 

•  44(4 

•  0094 

61 

.4574 

•  3(64 

lh 

.3.21 

,  .6512 

13 

.1416 

''  .3C1S 

3 

3. COCO 

.2414 

) 

4, 3>;C 

• - —  *5441-- 

3 

4.4563 

.6324 

9 

4.3iC0 

.0822 

; 

9.390ft 

.433) 

1 

3  .'.43  4 

»  VitV€ 

HO.  0 r  o*  »* 

PlUN$t! 

-.00  JO 

9.4C6C 

-.5377 

J 

e.C3CC 

-.417*. 

J 

1.5444 

-.0421 

2 

3.C303 

-.261* 

3 

1.6132 

-  .  *  **631£. 

6 

3.4638 

-.4017 

m 

.3576 

-.363% 

76 

•  A31 

-.4036 

<4 

.657% 

-.2664 

72 

.9*74 

-.0556 

76 

S.I8SC 

- .iili  . 

72 

.567V 

•  39C* 

-50 

.6679 

•  8365 

76 

*  .1471 

.0317 

4 

‘  .?:?% 

•  0616 

1 

1. 4444 

•  8912 

3 

,  0.3566 

- - -  -3*21.  _ 

4 

0.6035 

•  CC24 

) 

4*1551 

.3577 

3 

2.383C 

•  3  444 

3 

3.5364 

tic  sact^i 


tic  tst/re*.  c 


„  ,  iV  „  ^  A,<>  ^ 


TABLE  3.  CONTINUES 


C*r-A*IS  OBJECT  EXAKKC  CASE 

tHA«  fLA**€  aiSIAhCt 

•  16«3?83F02 

x  mt*cc  9. 

» 

AVECXO: 

-.199A9C8E-3J 

.t-2S7fe2E-C2 

9. 

•t9W6a€-83 

.lA4S7i2£-«2 

0. 

-.1666755E-83 

•  62t*>S92E*C3 

t» 

- - - 

-  ^47&tr9JE-)2 

3. 

.IAC9J\A*-CA 

-.59905456-83 

0. 

-».X6«9«*f -3* 

*.4994* A5E-23 

8. 

.19  3WC6I-33 

-.X232i<.gE-02 

9. 

-*i*34MU>83 

-.12VI-4C-32 

2. 

.29922S6C-CJ 

-.3-<S419E-C3 

6. 

■■♦AW?f4(‘4t  .. 

--*v±SA;3£>;3 

8. 

.J865694E-8J 

-.32337A«»£-C3 

8. 

*.J229/fci.£-C3 

8* 

.22162 396*33 

./:jr**2C-j«. 

5. 

-♦22U2J9E-04 

•  7l6M4?t-Ck 

9. 

.$*995m-C<. 
~*4u»u:t«ii. — 
.279i/69*-2*» 

-  ^rau  tat- ci.- 
.22167261 -W 
••22L6724C  -83 
-,fe*6S7fc|£-93 
^7«74U.)1. 


-.195354*e-n 

.J9S3$tfcE-43 

.«%eJ2ae-n. 

-.-.•dttLS23E- JW 

.fc93*222E-H5 


.65398756-33 
»*4'I  J375E-0  X 
,55-%52JC-C: 
-.JS4WS2U-51 
.7/982276-:* 
-  •J>Z$1221L~2  3 
.$*n»7t->4 
?.S4JlS.7f-3* 
.356I775E-3J 
-.35&i725trJ3 


-.  m.27j£-3*> 
-.:«.*5JAJt-33 

-.tKVSJi3C-53 

.2C1*/56C-C2 

.24147S6C-3? 

*1>929I9£-C2 

..-.X2U91M-C2 

.2J51A97E-8J 

.2itIV97C-ii 

-./M/J7CC-&J 
.•./4i?jr3E-33 
-,1S2JMJC-01 
.-•.4424-.JI  *£3-  -  , 
-.S?498*8£-C3 
-.62k51W«>33 
-.?2f<.63JC-tJ 
-.22tA613£-a3 
-.l*ft3225£-t! 
--a3C<3225C-CJ 
-.>7ii.9J3f-JJ 
-.?/«•-  9l3t-C3 

.  |ii7J$7t«$! 
.*4i/JS7£-?8 


c. 

3. 


7* 

*. 

8* 

C. 

8. 

3. 

C. 


0. 

3. 


-.943S865E-5J 

.35S.2557C-13 

3. 

~ - V3635463£i33  . 

.}*t2S97£~£? 

3*. 

-  .  . "  ,  . 

-.7996i*U-51 

.I97293t£-C« 

3, 

•73'Mltl£-33 

.urjgjxe-cj 

5. 

-.t72329^£-CJ 

-.♦9992  92E-CJ 

C. 

•l72J2J<.t-3  3 

».<&9&3«2£-33 

3. 

•:i9A966C*C3 

-,S7yi5lvC-3J 

e. 

- ~.2t44tb6£-JI_ 

-.627X91At-33 

3. 

,  v  .  M «,, 

.29<n640E'(J 

-.535639*6-53 

c. 

-^99:SAi£>4J 

r.56'.tlW£*Cl 

3. 

.I7SC9295-93 

-.<»23AJb6£-:i 

0. 

-.!255329t-CJ 

-.w23S3fc6£-33 

3. 

-.2«»>49i;E-3e. 

-.JA39997C-CJ 

9. 

.  4.UJA*37C’83 

-8* 

►Sl*30?*£-C? 

.599*356t-CJ 

.?«SlA2&9b£-tX 

.38V*J56£*CJ 

3* 

-.879J17jr-i3 

.♦*«.76Hf-;3 

K 

fcafer. - rtry ■  r  i &  ilrJ  VjVr'ili^^lir;yr,i ^  - 


TABLE  3.  CONTINUED 


.*Z*L3t9l-ZS 
-.i**4952C*«2 
.:n«e?f.}? 
-*7rtMK»  3J 
--  .^AM4E*ii 
-.5ll«20£-3J 
cr - .*444f**$-3*. 


~.3&3M3tf'»0t 

* - .J43iJ2iO#k 

.l?2J855?-«i7 
- »S?2.S*S$£-3? 

-.M?*929C-8* 

•»22077Ji£-G3 

.mmi£~43 

.?3J4*33£-JJ 

-.2anv*#ir-03 
.1J*9*9JC-C2 
- m»4MU»C? 


.  'J?Ck3iC-C4 
•  cl7t  .UC-CJ 
•»23WM-rt 
•w^43J2kE-oi 
«-.%*esj73£«ex 
•.44I44I$C*S3 
-.•9j;755r-c3 
-•4937755C-JJ 

*.i77!97i£-63 

.lC53736£-t? 

•ism  rm-i? 

«»*6m?*£-Ct- 

•.«6«74/4f-C* 

•l*11142C-S3 

-.l*:3l*2£.CJ 

-.l3bJJ?l£-:.« 

-.10US2IC-C? 

.wbl>635£-33 

•*4J/4*5C-4J 


<W  5I2£ 
.J343J31Jf-32 


fit  sucivj 

J3I33JI-;  J 


ttc  TIC  c«t«ctcts 


5»tfl!»S  AS-'Q» 


Cii.ti  r»ei» 


t*«£t  enitMitioi. 


tic  £a?hAS25« 


»IC  CHicjctts;. 


^•asi  iv:u‘ 


*  a. css 

1.3356 

5C.3C2 

•  9fe*5 

<44.53; 

.9487 

£56.33) 

234. *34 

25j.;>> 

.72*1 

JBC.iiO 

.4346 

tss.ccs 

.6^5 

*63.355 

•  *6«« 

*5  2*3)3 

.4  J  Jfc 

6K.dC 

.,1449 

556.239 

.4323 

4S0.4S3 

.>70* 

459.3SC 

•  2V« 

733.3)3 

.?J9? 

756. 33C 

.71*5 

f4c  :c; 

.  1 9*3 

456. SCC 

•  1591 

133,32) 

.£553 

-''■*  r—  *-a-«*  *-** 


TABLE  3*  CONTINUED 


*r*€9*  19 

*T6»T|I6  9*\9» 

0.3336 

TXC  E**>*fAStS« 

*  IIC  C*15*C 

:*>s*  - 

r*t<l**Cv 

•»?r 

pMtse  **. 

3.3)3 

1*3333 

0.9899 

56.610 

.*:?9 

-.9316 

toe. see 

.6371 

-.032? 

l  3*399 

.SJ<S 

-.0**7 

*3'»2J3 

.6*76 

-.1305 

?5*  363 

*5775 

-.219. 

333.322 

.-766 

-.37*7 

251 • SCO 

.<•677 

-.5726 

-36.923 

*U18 

-.7654 

45). 35) 

.12*? 

-.9*11 

. »(*I0S  -  . 

.2*75 

-1.3076 

555.323 

.267* 

•2.2205 

416rCC3 

.2*25 

-1.3503 

*  552.333 

.223 2 

-1.4/3) 

ns. 363 

•2C7J 

1.5694 

no. so? 

•  19f7 

1.5:36 

93C.CS3 

aic: 

1.4613 

056.899 

.17*7 

1.4953 

633. SCO 

.S6l7 

1.4953 

0 £lU  F®?Q*  55.)))) 


IttCET  TOT  AT  104* 


HC.  V  WIPSt 


i»*?tfU'C  a:ot*» 


•.3K> 


AP£*tU»£  >£I5Ht« 


*t»ixas*  -e.mo 


-•  r»u?u:c» 


-i.iJti 


T||,T  •*Ct5« 


*PCft*U3£  *<V* 


TABLE  3.  CONCLUDED 


s  YU«JC 

MC.  CT  ?i'S 

XNtCXSKfY 

-V»i 

3 

3.0030 

-.3027 

? 

3. aid 

5 

*.**&$ 

-.3621 

1 

3.33(3 

•.:»! 

1 

•  323) 

c. 

.mo 

-.oat? 

9 

.1133 

-.3991 

23 

•  tier 

-.ecofe 

n 

.74)3 

-.cos? 

•9MC 

-.CC!c 

50 

1.0CC5 

—  -  .1323 

64 

.*4:c 

•  933c 

.  !• 

.764  0 

♦  33)3 

23 

'6>« 

.331? 

9 

.11(3 

.35  XS 

<* 

•me 

•ecift 

1 

•  62CC 

♦lie;  — 

3 

3.332) 

.8C  26 

3 

M3K 

.30  27 

S 

3. SQIS 

•  833C 

3.33C4 

r  wlu€ 

SC,  C*  «T$ 

2«.T£*?;tr 

3 

3.0630 

».0*?7 

3 

3. 0*63 

-.3C?i 

3 

MJ1J 

-  -.»C21 

-  3 

>• )3)2 

-.)3ie 

2 

.C600 

-  -.1015 

.1200 

-.CC12 

"? 

Ate: 

-.3119 

10 

.7203 

-• 2)3f  * 

V# 

.V43C 

-.seal 

V? 

U3C93 

-.210; 

53 

:.S£K 

•C3*J 

56 

.760 

•  ))3t 

*6 

.4*51 

.5315 

t* 

.2*30 

.061? 

11 

.?3«1 

.0025 

5 

.1501 

•oeif 

•  S24J 

.**i\ 

<• 

.(191 

.55?6 

2 

.3633 

.3027 

2 

•  0633 

•0314 

5 

3.C35J 

tit  s .5>j?  tie  .••.■ujttrt.  «  e»cictt*s< 


m 


i 

1 


38 


^0.002- 

(<  ‘  '  .  . 

'«  - 

/\  ’  ,  s. 

»  „  — i 

^  - 

<■ 

*  / 

46.0W.C.  •  ■  it: 

-0,003 

;  r  "i  :•  V.  -- 

-0.002  -0.001 

j 

». 

«.  «  :  i  i  i  i  -  t  t  i 

0  '  0.301 

Z-AX1S 

i  i  Jo  -i 

0.002 

SLUR  CENTER 

X  «  0.00000  «  IC*1 

V  -  O.OOv'OO  *  10*01 

.< 

WAGE  PLANE 

2  -  1.50000  *  10*" 

Figure  20.  Y,Z  ray  profiLe  ior  on- axis  object. 


40 


•jit  taSLtaii- 


1.2. 


as< 


"  -  THETA(I) 


2 

,5  .  (, 

P 

< 

^  0.$  f- 

3  v 

o 

S 


* 


'  * 

A 

'&{  V5*% 
/ 

'  ~  Vi 


04. 


^  i. 


0t 

0 


200 


l 

<00 


600 


FREQUENCY 


600 


.0 

1000 


Figure  22,  Modulation  transfer  function  for 


on-axis  object. 


42 


v*V',<r_  .*■»  r  v  - 


C.002-  ' 


*  n« 
,j t  O'* 

n‘- 

>•  . 


■-{iV  if*  s * 


-o.oox- 


— 0.003L. .  /  . .  -i  j  .  /  .  .  (  ; 

-0.003  -0,002  -0.001 


-  r  . 

o 

X-AXIS 


*  •  4  *  *  •  t  I  . 

0.001  0.002 


BLUR  CENTER  X  -  0.00000  *  10'01  Y  «  2.62497  x  10-01 
IMAGE  PLANE  Z.*  LSOOOO  x  lO^1 

Figure  24,  Spot  diagram  for  off-axis  object. 


MODULATION 


1.2 


K 


i 


i 


1.0* 


A 


0.8  -  5  * 

•i  I 

-  h  ■  THETA(I) 

a  o  -  THETAI2) 

0.6  .. 


0*- 


0.2- 


0*. 

0 


200  400  ‘  600 

FREQUENCY 


i 

800 


>000 


Figure  27.  Modulation  transfer  function  for  off-axis  object. 


no  non  non  '»  o  on 


Appendix  A 


FORTRAN  LISTING  OF  MACE  ANALYSIS  PROGRAM 


PROGRAM  MAIN (INPUT, OUTPUT, FILKPL,TAPE5=I«PUT,TAPE&=3UTPUn 
COMHON/RAYS/X(20ll*Y(20i>,Zi2011,QX(2Qll,!m20t>,(m20i», 

1NPR, LABELS-?  ,XfiVG»rftVG,0£l.,2 
DATA  ICSNT/5HCENTER/ 

REAO  INPUT  PARAMETERS 
i  RE  A3  (5 , 13GHL  ABEL  <  I)  ,  1  =  1 ,5) 

IF (EOF (5) l 13, 11 
ID  STOP 
11  CONTINUE 
103  FORMAT (SAID) 

WRIT£<fe,lO?>  (LA3ELU)  ,1=1, »» 

1U7  FORMAT  (1M1,SA10) 

QELO  IS  SIGNF.O  DISTANCE  FROM  LAST  SURFACE  VERTEX  TO  IMAGE  PLAKF 
R£AO(S,102»OEL9 

102  FORMAT (F 10.01 
MP.ITE  (S,lil9) 

MRITE(6*13fl)0EL0 

103  FORMAT  (1H  ,  2GHIMAGE  PLANE  OIST  ANC.: » 5XE13. 71 

••**»«••••*  ****»•»«* 

RE AO  DECK  PUNCHED  SV  FOLOP 
READ (5,103) NPR 
103  FORMAT (I1G1 

00  10<*  1=1. NPR 

10*  RE AD (5 ,135) X ( I 1 ,Y (li  , 2(1) , On ( t! ,07(1) »GZ<I) 

105  FORMAT (a£l3.7) 


GOMPUTE  COORDINATES  ON  IMAGE  PLANE 
DO  iC6  1=1, NPR. 

icaisioxdJ/OiZUimioELn-zaiMxci 
y  ( i)  =  ( on  i » /or  u>  i  *  (  oeld-z  ( i  n  at  a  i 
ISs  Z(I)=M. 

FHPR=NPR 
XAVG=0 » 

YAVG=0. 

TRANSLATE  ORIGIN  TO  CENTER  OF  BLUR 
00  201  1=1, NPR 
XAVG=XA*S+X(I1 
201  YAV5=YAVG*Y(I) 

XAV3=XAVG/FNPR 

YAVG=YAVG/FMP« 

00  202  1=1, NPR 
X(I|  =X (11 -KAVG 
c 02  Y J I) =Y ( 1 1 -YAVG 
205  HRITE(o, 1031 
103  FORMAT (1H  I 

HRITE(o, lltf) XAV5, YAVG 

llO  FORMAT (1H  »3HX  AVERAGE ,5XE 1 3. 7 , 10X , 3HY  AVERAGE .5XE13, 7) 
NRITE(o, 1031 
00  203. K=1»NPR 

203  MRIT£(?m20.,)X(KI  ,YOO  ,Z(K) 

20*  FORMAT (1H  ,3(5XE13.71) 

CALL  INITIL 
CALL  RAPLOT 
CALL  MTF 
CALL  SCAMS 
GO  TO  1 
cNO 
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suoroutine  raplot 

Ci3fHGN/RAYS/X{201)  , Y C 201) .2 <20 1)  ,QX  1201) .  QY  {201)  *  dZ<  20t»  . 
lNPR,LA3ELi5)',XAVG,YAVG,rELO 

R£A3  IfPUr  3ATA 

READ  <5^  1G3).SCAL£»  OELT  XY,  NENPHS  ,NCHRS 
153  FORlAT<2FiO.O,2IlC) 

WRIT£<o,t0v> 

I0i»  FORMAT  UP1.75H  PLOT  SI2E  TIC  SPACING  TIC  EHPH 

1ASI2EO  rIC  CHARACTERS) 

WR7T£< 6,105)  SCALE, DEL  TXY,N£><?HS,NCHfiS 
105  FORMAT  { iri  ,2<4XE15.d)  ,2U5XIh)  ) 

XLEFT=-SCAl£ 

XRIGHT=SCALE 

OSAW  X-Y  P..3T 

CALL  GRIOULEFI  «XRIGWT,Xt  EFT  BRIGHT  ,OELTXY  .OELTXY  ,NEHPHS,  NEKPHS  ♦ 
1N£MPHS,N£HPHS,«CHRS,NCKR3) 

CALL  PRINT/<60,LA8£L, 270, 1000) 

CALL  PRINT  Yl -6  *6HX-AXIS»  530  *  3) 

CALL  APRNW<L,-lf«t-S.  6HY  AXIS, 0,550) 

CALL  APL OT VIUPR, X,Y, 1,1, 1,44, IERR) 

CALL  PRINT  Y( -37 ,37H3!.'JR  CENTER  X  a  V  =,100,20) 

CALI  LASLVUAVG,2-»4»2Q,-6,1.1> 

CALl  L43LVtYAVG,404,20,-6,l,l) 

CALL  PRINT V<-1 7,1 7HIH AGE  PLANE  2  =,1C0,0) 

CALL  LA3LV(0£L0, 244,0 t-6,1,1) 

OR AW  Y-2  Pl3T 

CAll  6»ID< XL£FT,XRIGNT,XLEFr,XRIGHT, OELTXY, DEL TXY,NE<JPWS,NENPHS, 
lN£M?HS,«£NPHS,NCrtRS-:OC)«S) 

CALL  PRINT  Y<60  »LA8£L  *  2? /J  *  i 0  30) 

CALL  PSINTY<-6,6H2-AXIS, 500,0) 

CALL  APRNT Y(3  » -!<•„  -6, 6HY  AXIS, 0,600) 

CALL  PRINTv'(-3?,37h8LUR  CENTER  X  =  V  =,100,20) 

CALL  LA0LV!XA'YG,244,20,-6,1,1) 

CALL  LA3LV(YAVG,4  0.4,20,-6,1,1) 

•CALL  PRINTY<-l7,17HIHA3S  PLANE  Z  =,100,6) 

CALL  LA8LVOELO, 244, 0  ,-6, 1,1) 

00  250  1=1, NPR 
I21=NX V{XLccT ) 

IZ2=NXV( KRIJHT) 

YL  =  Y  1 1 )  +  XLEcT*QY ( 11/02(1) 

.  IF(A8S(Yt) ,l£ • ASS  CXLcFT) )G0  TO  201 

ZLSS«=Xi.Eff*<SIGW<XLEf*,YU-ym)/(YL-YlII  ) 

YL=SIGH(XL£FT,YU 
I?C=NXV{ ILnEW) 

201  YR=Y <I‘*XRI3HT*QY<I>/QZ<I) 

IF<A3S(YR).LE.A8S(XLEFT>)C3  TO  202 
ZRNEW=XRI&MT*(SrGN(XR IGHT »  YR)  -Y<  I)  )  /  {YR- Y  ( I  > ) 

YR=SIGN{XRI3HT,YR) 

I22=«XV(2RNEH> 

202  IY1=NYV(YL) 

IY2=i)YV(  YR> 

200  CALL  H!*£VUZl,iYl,XZ2,tY2> 
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OR  AH  X~Z  Pl.1T  o 

CALL  GRI0(XL£FT»XRI3HT«XL£Fr » XRIOH  T»Q£LTX Y  »OELTXf  *N£8PHS»l«£HPHSt 
INEHPHS*NEH*MS,NCHRS»NCHRS> 

CALL  PRXNTV(6G,LA3EL, 270,1000) 

CALL  PRlNTVi-{f,SHZ-AXIS, 530,0) 

CALL  APRNTtf (0 ,  ~1<,  6HX  AXIS, 0,600) 

CALL  PRlNrv(-??,C7H3LUR  CENTER 
CALL  LAaLV(XAVtC,2**4, 20,-6, 1,1) 

CALL  LA3LV<YAV&,<,04,20,-6,1,1) 

CALL  PRI NI V(-l?,l 7HIKAGE  PLANE 
CALL  L  A^iLV  OELO,2 LA,  0  ,-3,1,1) 

.00  300  I=1»NPR 
I2l=NXV<XL£PT> 

IZ2=NXO<XRI3HT) 

XL=X<I)+XLt“T*QX(I)/3Z<I) 

IF(ABS<XL)«L£, ASS (XLEFT) )  50  TO  301 
ZLN£H=XL£FI*  (S  ION  < XLEFT  »XL)«X(I>)/  (XL-X ( I  > ) 

£L=SIGMXLcFT*XL) 

IZl=NX Y ( ZLNEX) 

301  xR=xa)»xRi;Hr*axa>/Qza>  . 

IFIABSCXRK^LE.ASSIXLEFT) >G0  TO  302 
ZRN£W=XRI3NT* (SIGN (XRICHT , X?) -XII)) / (XR-X C  I)  > 

XR=SIGN(XSI3HT,XR> 

IZ2=NXV(  ZRNZW)  i 

302  IX1*NYV(XL) 

IX2=NYI/(XR> 

300  CALL  LlN£VUZl,IXl,IZ2.IX2l 
RETURN 
ENR 


x  *  r  =»ioo,2e) 

Z  =,100,0) 


ooooo  ucoo  <..  o  o  u  uuu  o  o  o 


SUBROUTINE  RTF 


A  SUBROUTINE  DESIGNED  TC  COMPUTE  THE  HOOULATIOA.  TRANSFER  FUNCTION  Of  AH 
OPTICAL  SYSTEH  WITH  MODERATE  ABERRATIONS. ALSO  COMPUTED  IS  THE  PHASE 
ANGLE  PHI  ASSOCIATES  WITH  HE  KTF. 

COHHON/RAYS/X 1 201 ) » Y ( 20l» , Z  <CS1 ) ,QX (201 » , <5Y (20 1 * . 3Z( 2011 .MUM, 

1LABEH6I 

DIMENSION  X?{5u)«VP(5&5vA>'rtrF<50)  ,F«£Q<50)  ,PHI*50>  ,TX(201> 

DIMENSION  METAC2) 

data  cvTR,MOPi/o.omsi£<),s.2e3i353/ 

READ  in  PLOT  CONTROLS.  TO.  Dr  FREQUENCIES,  STARTING  FREQ., 

OELT A  FREQ. ,  AND  TWO  AXIS  ROTATION  ANCLES 

READ  15,110)  NENPrtS.HCHR!  .NFREO.STFKEQ. OFREQ, THETA (1) .THETAI2) 

10  FORMAT  f-3  II  J,<.Fl0.  4) 

TNUM=NUh 

SET  UP  PLOTTER 
A1=STFREQ 

A3=Ai,(rLOAT(NFR£Q>*l.t<»*OFREQ 

A2=0. 

A*=l.? 

CALL  GRID t Ai, A3,A2,A<, ,CFRE3« . 05, NYMPHS,*, , N£NPHS,<f ,NDHRS, 3) 

CALL  PSINTi/f&3, LABEL. 2.  0,1000) 

CALL  P3INTYI-R,3MFREQU?HCY, 500,3) 

CALL  APRHTVt  0,  -14 ,-10 ,  ‘-QHMODUL ATI0N.3, 600) 


DO  2J5  L‘»=l,2 

FREQii)=STFREQ 

ROTATE  AXES  USING  THETAtLS.)  AND  COMPUTE  f'PRlMEiKt  rOR  EACH  X(K). 
OO  115  <=1,NUM 

5  TX<O*X<>i),C0S<TH£TACL‘.)*CVrR)  ♦Y(X)*SlNtTHETA(L.,l  *CYTR) 

00  130  J3=1,MFRE  0 
REG1-G ,3 
R£G2=0.0 

>  00  175  J2=1,«UM 

REOi =R£G1*C3S ( TMQPI*FR£Ql J3> *TX ( J2 ) ) 

REG2  =  RtG2fSINlTHOPI»FRE5;U3)*TX(J2>) 

135  CONTINUE 

AS=REG2/TNUM 
AC=RtGi/TNUN 

COHPUTc  MTF  FROM  AS.AC.ANO  MUM. 

ANHTF (J31 =S1RT{AC*AC»-S*AS) 

PHI { J3) =ATAN <  AS/AC) 

XP<J3)=FRSU(U3) 

YP<J3) =ANMTF< J3J 
FR£OIJ3H)=FREO(J3)«OFREQ 
Uu  CONTINUE 


52 


c 

c 

c 


PRINT  OUT  HARO  COPY  OF  RESULTS. 

HRIT£(o,lSj>  NFRE  :,STFREQ,0?R£Q, THETA(L4> 

150  F0RMAT(1H1,3X,*  NFREQs  *,I5,5X,*  STARTING  FREQ=  *,F10.4,6X,  * 

lOELT'A  FR£Q=  *,F10.4,6X,*  TARGET  ORIENTATIONS  *,F10.4,///> 
HRITEta, 160)  NEHPHS»NCHRS 

160  FORMAT (' H,3X»*  TIC  EMPHASIS*  *»I5,8X,*  TIC  CHARACTERS*  *«I5,///> 
HRITEJo, 155> 

155  FORMAT (10X,*  FREQUENCY  *,10X,»  MTF  *,10X»*  PHASE  AN3LE  *,/✓? 

WRIT£<6» 15$) .  <FREQ(J4>,ANMTF(J4i ,PHI(J4> » J4=l .NFRE3) 

150  FOR«AT<8X,Ft0.3,6X,FA0.4,8X,F10,4> 

C 

C  PLOT  MTF  VS.  FREQ. 

C 

IF<L4.EQ.1>  GO  TO  102 
IF (L4»  cQ. 2)  GO  TO  103 

102  CALL  APLOTV(NFR£Q,XP(ii,Ya<l>, 1,1,1, 63, IERR1) 

'50  TO  10  4 

103  CALLAf'LOTVlNFREQ,XP'i) ,YPti> ,1,1, 1,38, IERR2> 

104  NFR=NFR£Q-1 

DO  101  11=1, NFR 

-  IX1=NXV (XP (II )  ) 

IX2=NXV(XPUI  +  1>> 

IY1=NYV( YP(II) > 

IY2=NYV(YP(II+1M 
101  CALL  LINEV(IX1,IY1,IX2,IY2> 

215  CONTINUE 
Rf  TURN 
ENO 


! 

i 

t 

< 

\ 


53 


OOOOOUUOOUOUO  OOOOOOOOOCJOOOO  0,0  o 


SUBROUTINE  SCANS 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCDCCCCC'*'  CCCCOi  sccccccccccc 
0  THIS  SUBROUTINE  IS  DESIGNED  TO  SCAN  ACROSS  THE  SPOT  0  V  'AM  ON  THE  IMAGE  C 
C  PLANE  OF  THE  OPTICAL  SYSTEM, USING  EITHER  A  RECTANGULAR  CIRCULAR  A PER •  C 
C  TUBE. THE  USER  HAS  THE  OPTIONS  OP  CHOOSING  THE  CENTER  OF  SC^N  COLROINFtTES  C 
C  HIHSELFtOR  USING  THE  DIAGRAM'  CENTROID  AS  SCAN  CENTER. IN  ADC.-TIOKVXE  MA“  C. 

S  EXAMINE  THE  DIAGRAM  ^.ROH  A  SPECIFIED  TILT  ANGLE  AND/OR  A  SKFiFXED  TRAN-  C 
5  SLAT  ION  OF  THE  COORDINATE  AXES. THE  OUTPUTS  ARE  THE  LOCATION  Of-  HE  CENTERC 
3  OF  THE  SCAN  APERTURE, THE  NUMBER  OF  SPOTS  ENCOMPASSED  BY  THE  APE '  JRE  AT  C 
3  THAT  POSITION, ANO  THE  NORMALIZED  VALUE  OF  SPOTS/POStTION  RELATIVE  TO  THE  C, 
3  MAXIMUM  NUMBER  OF  SPOTS.SCAMS  ARE  PERFORMED  IN  THO  MUTUALLY  PERPENOIC-  C 
C  ULAR  DIRECTIONS.  C 

CCCCCCCCCCCCGOCCCCCCCCCCCCCCCCCCCCCCCCCCCOCCCCCCCCCCCGCCCCCCCCCCCCCCCCCCCCCC 

COMHON/RAYS/XI241) ,V(2fll) ,Zt201) , QX (20 1) ,  QY (20i> , 0Z(  201  > . NUM. 

1LA3EL(6> 

DIMENSION  YC<100>,,NNUMC<100>  ,  ANNUMC (10 0 ) ,  VNORM (108) 

OIMENSIQN  ANUHC (100) , ANORMt 100 ) ,XC(1U0),NUHC(1QQ) 

OAT*  REVOLV/6HROT  ATE/ 

DATA  ICIR/oMCIRCLE/ 

OATA  CVTR/. 31745329/ 

ANUH=NUM 

REAO  IN  APPROPRIATE  PARAMETERS  AS  POLLOMS— 

NS-NUHBEH  OF  STEPS  USEO  IN  SCANNING  PROCESS. 

SS2=VERTICAt  DIMENSION  OF  RECTANGULAR  APERTURE. 

SSI* HORIZONTAL  OIHENSIOM  OF  RECTANGULAR  APERTURE. 

R=IF  SPECIFIED, LET  SS1=SS2=0. RADIUS  GF  CIRCULAR  APERTURE. 

VTR*IF  CENTROID  NOT  USEO, GIVES  Y-COORQINATE  OF  SCAN  CENTER. 

XTR= IF  CENTROID  NOT  USEO, GIVES  X-COORDINATE  OF  SCAN  CENTER. 

TILT=R0TATE  COORDINATE  AXES  TO  THIS  ANGLE  ANO  COMMENCE  SCAN. 

K£Y=TYPE—  CIRCLE  —  IF  CIRCULAR  APERTURE  USED. OTHERWISE  BLANK. 
A1=NEGATIVE  VALUE  OF  SCAN  LIMIT. 

A2=P0SITIVE  VALUE  OF  SCAN  LIMIT. 

ORIENTsIF— ROTATE— TYPED  ON  SEPARATE  CARl>,$Sl  =  S32  ANO  SS2=SS1  FOR 
Y  AXIS  SCANNING. 

QX=TIC  SPACING  TO  BE  USEO  ON  PLOTTER 

ON*L ABEL  and  EMPHASIZE  EACH  NTH  TIC  ON  PLOTTER  OUTPUT, 

R£AD<5,700)  NS,SS1,SS2,R,YTR,XTR,TXLT,KEY 
TOO  FORMAT (I10,SF10.0,A6) 

MRITE(6, 79o>  NS.SS1.SS2, R 

796  FORMAT (lHi»*  NO.  OF  STEPS*  *,U8  ,5X,»  APERTURE  MIOTH*  *,F10.4, 

15X,*  APERTURE  HEIGHT*  *,F10.4,5X,*  RADIUS*  *,Fl0.4,///» 

HRITE(6,798»  YTR ,XTR , TILT , <E Y 

798  FORMAT (1H0 Y  TRANSLATE*  •,F10.4,5X,*  X  TRANSLATE*  *,F10.4,5X,*  T 
1ILT  ANGLE*  >,F1Q. 4,5X ,<■  APERTURE  KEY*  *,A6,///> 

11  REAO (5,701)  A1.A2 
'  701  FORMAT (2Fli). 0 ) 

WRITE ( 6, 802)  A1.A2 

„  802  FORMAT (1H0»*  LEFT  SCAN*  * ,F10 . 4, 10X,*  RIGHT  SCAN*  *,F10.4»///> 

RSA0(5,815)  ORIENT 
815  FORMAT (A6) 

IF (ORIENT. EQ.REVOLV)  WRITE (5,6231 

820  FORMAT (1H0»*  SCAN  RECTANGLE  IS  ROTATED  THROUGH  90  OEG  FOR  Y  SCAN 


oooo  0000000  \  o  o  o 


ii SteswSa 


-=.  W^’Wm 


-■vw  ywA'^A<‘'"^'\'  *  4v».-e 


820  CTO 


IF{TILT.Nc.3.)  GO  TO  3 

IFUYTR.EQ.O.) .ANO. (XTR.eQ.3.) ‘  GO  TO  10 


PERFORM  APPROPRIATE  COORDINATE  TRANSFORMATIONS  USING  ROTATION  ANO/OR  TRANS¬ 
LATION  PARAMETER.  -  ' 

i  TILT=TILT*CVTR 
SI1=SIN(TIlT) 

CIT-COS(TILT) 

00  15  J=1 *NUM 

Tl=< (J) -XT  R 
T2  =  Y ( J ) - YTR 
X(J) =Tl*CIT*T2*SIT 
L5  Y(J) =-Tl*SIT+T2*ClT 

L0  IF  (KEY  .EQ.  ICIR)  SS.t=SS2=2 . 0*R  CIRTEST 

HORIZONTAL  SCAN  BEGINS  HERE. 

COMMUTE  STEP  SIZE  FOR  SCANNING. 

ANS=HS-1 

SS= { A2-A 1) /ANS 

CO.* NT  NUMBER  OF  RAYS  >  HITHIN  THE  APERTURE  AFTER  EACH  SCAN  STEP. 


GO  120  J=i,NS 
3J= J 

XC(J)*A1-MRJ-1.0)  *SS 
NUMC ( J ) =3 
□  O  lie  <sl»  NUM 

IF ( KEY  . EQ.  ICIR>  GO  TO  112  CIRTEST 

XBL*XC (Jl -SSl/C.O 

XBR-XC ( J ) +351/2.0 

YBU=  +332/2.0 

Y30=  -S32/2.0 

IF(X(K).G£.XBL.ANO.X(IO.LE.X3R.ANG.Y(X>»5E»Y0D.ANO.Y(K> .CE.YSU) 

1NUMC ( J I =  NUM3 ( J ) *1 

GO  TO  US  CIRTEST 

112  .  C=SQRT  I CX  f<l-XC{J))**2t<Y  (Kl')**2>  CIRTEST 

'  IFJABS(O)  ,LS.R>  NUMC  T-JT^NUMC  (J>  +1  CIRTEST 

118  CONTINUE 
120  CONI INUE 

FINO  MAXIMUM  VALUE  OF  NO.  OF  RAYS  ANO  NORMALIZE  EACH  NUMBER  TO  THIS. 

MVAL=NUMC<i» 

30  13p  J=2.NS 

130  MVAL*  MAZO (MV At .NUMC <  J) ) 

00  .140  J2  =  l *  NS 
ANUMC  <  JR )  =N‘JHC  <  J2  J 
AMVAL-MVAL 

l*.0  ANORMt  J2>  *  ANUHC  (U2 ) /A  »Al  ' 

MRITElo. 203) 

200  FORMAT (IH1»3X»*  X  VALUE  *,13X.*  NO.  OF  RAYS  *S10X«*  INTENSITY  *♦ 

1///) 
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o  U  CO  'J  ooo  ooo  ooo 


NRITE<t>,?L:>>  (XC(J)rNUHC(U)«ANORM(J)  ,J=i,NS> 

205  FOR1Ar(6X,FiC.4,5X,UGUlX,-10“» 

VERTICAL.  SCAN  BEGINS  HERE. 

NOTATE  SCANNING  RECTANGLE  THROUGH  90  DEG  If  THIS  0PTI0H  HAS  CHOSEN. 

IF  (ORIENT  •  E'J.REVOLV)  GO  TO  207 
GO  TO-  208 
g 7  SS2?=3S2 

SS2=SSi 
SSH--SS2P 

COUNT  N'UMBER  OF  RAYS  HITHIN  THE  SCAN  APERTURE  AFTER  EACH  SCAN  STEP. 


>08  .00  220  J=i»NS  • 

:<U.=  J 

YC(J)  =  Ai<-,(RJ-i.O>  *SS 
NNUHC(J)=0 
OO  218  K*i.NUH 
,IF(XEY  .EQ.  ICIR)"  GO  TO  212 
X8R=  +SS1/2.0 

:X3L  =  -S31/2.Q 

Y8U*YC(U)  +S32/2.0 
YBO=YC(J1'-SS2/2.0  : 

IF<X(K).GE.X3L.AN0.X(K).LE.<3R.AN0. Y(K) .S£ ,TBO.ANO.Y(K> .LE.YBU) 
INNUHC(J) =NNUHC(JJ+1 
GO  TO  218 

212  QrSURT  ( ( Y  <K)-YC{J))**2*X  (0**2) 

IFtASG(O)  .LE.  R)  NNUMC( J) *NNUNC ( J) +1 
218  CONTINUE 
22u  CONTINUE 


CIRTEST 

-CIRTEST 

CIRTEST 


FIND  MAXIMUM  VALUE  OF  NO.  OF  RAYS  AND  NORMALIZE  EACH  NJHBER  TO  irtlS. 


KVAL  =Nr,'UHC(ll 
00  230  -J=2.NS 

HVAt=HAXO(ftVAL.NNUMC((J)) 

00  240  J=1,NS 
ANNUNC ( J) =UNUNC( J) 

„  AHVAL^HVAL 
2«.0  VNORH(  J)  sANNUNCl  J1  /AHVAL 
HRITE(G,30j) 

300  FORHAT (IH1*3X»*  Y  VALUE  *,10X»*  NO.  OF  RA>S  *,10X,*  INTENSITY  *f 

\//n 

HRITE(o,3CG)  (YC(L) «NNUMC(L) .VNORN(L) *L=lt>'S  ) 

305  FOR.1AT(9X,F10.4,5X,I10.11X,F10.4) 

PLOT  RESULTS 


A3=0. 

A4=1.2  ,  , 

REAO(S,430T  OX.N, NCHRS 
400  FORMAT (F10.4.2I10) 

WRIT  E (6»  8O4)  OX,N,NCMRS 
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...  .  i-.- : y.V **¥ ■>  '  ■■ 

',  44 


FORMA!  U81,*  TIC  5fACING=  *,fiO«4,BX,»  TIC  EHPKASIZiD*  *,Ifri<U*. 
1*  Crt ARACT£RS=  *,I6: 

CALL.  GRIOUi,A2,A3,A4»CX..03,Nf4fN»4,NCHRS<3) 

CALL  PRINTVibt) , LABEL ,270, 1030 > 

CALL  PRINTX<~9,9HSCAN  AXIS, 500., 0>. 

CALL  APRNTVtt!,“14.;-9,9HINT£9SlTY, 6,8003 
CALL  APLOW<NS*XC<i>*ANda>mftltV*.lfS5*ItRR.3) 

MSS=NS-i 
OO  410  JJ=i«NSS 
JXi»HX¥C«C-Uj)  > 

IX23NX V  C  XC.L  jtJ+i) ) 
m*NYV(«N02M<JJ> ) 

rr2=«Y¥(ANJRH(  JJ+DI 

CALL  LIHEV<rxi,ITi»£X2.IY2! 

CALL  APLOTtf iNS  ,YC!  1)  .VNORMUI  *  1*1 , t ,  38c  TERRA  i 
NSPP=NS  -1 
00  420  >0<=1..JSP? 

IXi=NXV(YC{<K» )  , 

ix2»mxv  fycc<K*in  '  . 

IYl^NYVC VNORMtKO  !  ,  . 

IY2=NYV(VN6l?H<J«+i)>  v 

CALL  LINeVUXt,mPiX2.,IY2> 

RETURN 

£N0 


SUBROUTINE  3KIO(XLfcPT  ,  XRiGHT,  VOOWN, YUPfOX, OY»K,M» I,J*NX,MY) 

ORAR  ANO  LA3EL  GRID 
CAIL  FftANE/«5»  c' 

CALL  XSCAlyTXLEFT, XRIGHT, 86,43! 

CALL  YSCALV(YG0HN,YUfi,89,43> 

CALL  LT.N£VC9«,.88,88,96£i) 

CALL  LIN£V<88,980,98fl,980! 

CALL  L  INEV,.(_980 ,  9SG  ,980  »,86J 
CALL  L IN£V 1 98, 88, 980 , 88) 

CALL  LINRYil *76,68, 98* XLEFT * XRIGHT ,0X*N,I , NX,3) 

CALL  lINRV (2,24, 68*98 ,  YOOWN, YUP«OY , M,J,NY,10> 

CALL  LINRV(i,?jD,9?a.98C,XLc"T,XRrGHT,flX,N,0,NX,8)  . 

CALL  L INRtf { 2  » 24,970 ,980 ,YOOWN» YUP*OY»N»3*NY,tB! 

RETURN 

END 


SUBROUTINE  INI  TIL  'Y 

INITIALIZE  PLOTS 
CALL  CA«RAY(9,9» 

CALL  FRAHEVTO) 

CALL  XAXISV10,0> 

CALL  X AX ISV (C , 1023) 

CALL  VAX ISif (G , 0} 

CALL  YAXISUd  023,  0 ! 

CALL  PRI«TYJ-2i,2IHPLOTS  FOR  HYNAN  CURRY, 360,490)' 

RETURN 

END 
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Appendix  B 


Format 


Card  Group 


7  (1  card) 


(110,  6F10.0,  A6) 


8  (1  card) 


(2F10.0) 


Parameters 

DFREQ  =  Frequency  increment 
THETA (1)  =  First  scan  angle 
(degrees) 

THETA (2)  =  Second  scan  angle 
(degrees) 

NS,  SSI,  SS2,  R,  YTR,  XTR, 
TILT,  KEY 

NS  =  Number  of  scan  steps  or 
positions 

551  =  X  aperture  dimension 

552  =  Y  aperture  dimension 
R  =  Radius  of  aperture 

YTR  -  Y  translate  coordinate 
XTR  =  X  translate  coordinate 
TILT  =  Axis  rotation  (degrees) 
KEY  =  "CIRCLE-'  if  circular 
aperture 

KEY  =  Blank  if  rectangular 
aperture  ~ 

Al,  A2 


* 

1 


9  (1  card) 


(A6) 


Al  -  X  starting  point 
(smallest) 

A2  =  X  ending  point 
(largest) 

ORIENT 


3 


10  (l  card) 


(F10.0*  2110) 


ORIENT  *  "ROTATE"  to  rotate 
Y  aperture 

ORIENT  *  Blank  for  no  rota¬ 
tion  of  Y  aperture 

DX,  N,  NCHRS 


OX  -  Tic  increment 
N  -  Nth  tic  to  be  emphasized 
and  labeled 

NCHRS  ■=  Tic  label  characters 
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